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1.0  INTRODUCTION 


In  the  past  few  years,  the  National  Highway  Traffic  Safety 
Administration  (NHTSA)  has  displayed  a great  deal  of  interest  in 
determining  the  underlying  causes  of  vehicle  rollover  accidents 
and  developing  vehicle  and  roadway  design  criteria  which  will 
reduce  both  the  number  and  severity  of  rollover  accidents.  Vehic- 
le rollover  is  expected  to  be  a growing  problem  as  NHTSA  1986 
accident  data  files  show  that  utility  vehicles  and  small  pickup 
trucks  have  the  greatest  number  of  rollover  deaths  per  million 
vehicles  (155.7  and  127.7  respectively),  while  also  having  the 
highest  percentage  of  fatal  accidents  involving  rollovers  (66% 
for  utility  vehicles  and  45%  for  small  pickups) . At  the  same 
time,  sales  of  these  types  of  vehicles  have  been  growing  dramati- 
cally with  the  market  share  of  utility  vehicles  increasing  from 
2%  to  6%  in  the  past  decade. 

There  are  two  types  of  vehicle  rollovers  which  occur 
frequently  in  single  vehicle  accidents.  A tripped  vehicle  rol- 
lover can  result  if  the  vehicle  strikes  a curb  or  other  obstacle. 
An  untripped  vehicle  rollover  occurs  during  an  abrupt  vehicle 
maneuver  without  the  help  of  a tripping  mechanism.  Both  tripped 
and  untripped  vehicle  rollovers  are  equally  dangerous  and  de- 
tailed investigations  of  their  behavior  are  needed  in  order  to 
determine  their  causes  and  thereby  reduce  vehicle  rollover 
propensity. 
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Several  projects  have  been  initiated  to  study  the  factors 
involved  in  rollover  accidents.  Calspan  Corporation  studied  the 
rollover  potential  of  vehicles  on  embankments,  sideslopes,  and 
other  roadside  features  on  behalf  of  the  Federal  Highway  Adminis- 
tration [1].  The  Tripped  Rollover  Vehicle  Simulation  [2]  deve- 
loped by  Systems  Technology,  Inc.  (STI)  is  being  used  by  NHTSA  to 
predict  vehicle  rollover  response  during  maneuvers  involving 
lateral  impact  with  a rigid  curb.  The  state  of  Maryland  is 
currently  collecting  accident  data  which  will  be  used  to 
determine  the  statistical  correlation  of  factors  associated  with 
vehicle  rollover,  and  the  University  of  Missouri-Columbia  is 
currently  developing  an  Advanced  Vehicle  Rollover  Model  to  study 
the  dynamic  behavior  of  vehicles  in  both  tripped  and  untripped 
vehicle  rollovers. 

The  STI  Tripped  Rollover  Vehicle  Simulation  can  be  used  to 
predict  a vehicle's  response  during  a maneuver  in  which  the 
vehicle  slides  laterally  into  a rigid  curb  with  some  small  head- 
ing angle.  The  STI  simulation  allows  the  user  to  select  dif- 
ferent sets  of  initial  conditions  such  as  lateral  vehicle  speed, 
angle  of  impact  with  the  curb,  initial  distance  from  the  curb, 
ect.  The  model  has  eight  degrees  of  freedom  which  include  vert- 
ical and  horizontal  translations,  roll,  pitch,  and  yaw.  While 
this  model  generates  vehicle  response  in  tripped  rollover  situa- 
tions, the  influence  of  various  parameters  on  rollover  propensity 
is  not  immediately  evident. 
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Knowledge  of  how  the  different  characteristics  describing  a 
vehicle's  dynamic  systems  and  sub-systems  affect  its  behavior  in 
accident  situations  is  an  essential  aspect  of  understanding  the 
rollover  problem.  However,  based  on  simulation  runs,  it  is  not 
always  clear  which  characteristics  have  the  greatest  influence  on 
the  behavior  of  a complex  mechanical  system  such  as  a vehicle. 
Also,  attributes  which  may  be  important  in  one  situation  may  be 
unimportant  under  different  conditions,  making  the  entire  analy- 
sis even  more  complex.  It  has  been  shown  [3-8]  that  sensitivity 
methods  are  very  useful  in  determining  the  influence  of  various 
design,  operational,  and  environmental  parameters  on  the  response 
of  complex  mechanical  systems.  As  a result,  it  was  proposed  that 
the  sensitivity  algorithms  developed  by  the  University  of 
Missouri-Columbia  be  applied  to  STI ' s Tripped  Rollover  Vehicle 
Simulation  to  gain  insight  into  the  influence  of  certain 
parameters  on  tripped  vehicle  rollover  behavior. 
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2 . 0 PROJECT  OBJECTIVES 


The  primary  objective  of  this  project  was  to  conduct  a 
sensitivity  analysis  of  STI ' s Tripped  Rollover  Simulation.  The 
sensitivity  of  vehicle  system  parameters  were  investigated  and 
ranked  in  order  of  importance  in  tripped  rollover  situations. 

Relevant  system  response  variables  such  as  sprung  and  un- 
sprung mass  roll  angles,  roll  velocities,  lateral  velocities, 
ect.  were  analyzed  using  percentage  sensitivity  functions.  The 
percentage  sensitivity  function  was  chosen  to  provide  the  clear- 
est physical  interpretation  of  results.  In  addition  to  system 
response  variables,  an  additional  function,  vehicle  rollover 
prevention  energy  reserve  (RPER) , was  derived  based  on  energy 
considerations  and  was  used  as  an  indication  of  vehicle  rollover 
stability.  The  sensitivity  of  this  function  to  system  parameter 
changes  was  then  conducted  using  various  sets  of  initial 
conditions . 
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3 . 0 


STI  TRIPPED  ROLLOVER  MODEL 


The  STI  tripped  rollover  model  [2]  simulates  the  motion  of  a 
vehicle  which  skids  laterally  on  pavement  and  comes  into  contact 
with  a roadside  curb.  The  model  simulates  the  vehicle's  skidding 
motion  prior  to  curb  impact,  as  well  as  vehicle  motion  after 
impact.  Vehicle  rollover  is  assumed  to  take  place  when  the  roll 
angle  of  the  unsprung  mass  exceeds  approximately  85  degrees. 

The  simulation  uses  a non-linear  eight  degree  of  freedom 
model  which  utilizes  simple  linear  subsystems  to  model  tires, 
suspension  and  impact  force.  The  vehicle  model  (Fig.  1,  3 and  4) 
consists  of  two  masses,  a sprung  mass  representing  the  vehicle 
chassis  and  a single  unsprung  mass  which  represents  the  combined 
front  and  rear  suspension  systems.  The  sprung  and  unsprung 
masses  have  separate  degrees  of  freedom  for  roll,  lateral  and 
vertical  translation,  while  vehicle  yaw  and  pitch  are  analyzed 
using  a single  lumped  mass.  The  impact  force  for  lateral 
wheel/curb  impact  is  modeled  using  both  plastic  and  elastic 
deformations.  Damping  effects  are  included  through  energy  dis- 
sipative forces  in  the  vehicle's  tires,  the  lateral  bushing 
between  the  sprung  and  unsprung  masses,  the  shock  absorbers  in 
the  suspension  system,  and  the  wheel/curb  impact  force. 
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Fig.  1 
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3 . 1 Tire  Model 


The  tire  normal  reactions  (^ZLF'  ^ZRF'  ^ZLR'  ^ZRR^  a.re 
based  on  a linear  tire  deformation  model.  The  tire  normal  forces 
are  proportional  to  tire  deflections  caused  by  vehicle  heave, 
pitch  and  roll.  The  tire  deflections  are  determined  by  the 
finding  the  difference  between  the  undeflected  tire  radius  and 
the  distance  between  the  axle  and  the  road  surface  under  each 
wheel.  Tire  damping  forces  are  also  included  and  are  proportional 
to  the  tire  rates  of  deformation. 

Tire  lateral  sliding  forces  (FyLpr  ^YRF'  ^YLR'  ^YRR^ 
functions  of  tire  normal  forces,  tire/surface  frictional 
properties,  and  sliding  velocity. 


3 . 2 Suspension  Model 


The  vehicle  suspension  model  consists  of  two  springs  which 
have  upper  and  lower  bump  stops.  The  springs  deflect  as  the  roll 
angles  of  the  sprung  and  unsprung  masses  change,  and  as  the  two 
masses  heave  in  relation  to  each  other.  The  bump  stops  prevent 
excessive  roll  angle  and  heave  differences  between  the  two 
masses.  Damping  in  the  suspension  is  included  to  represent  both 
shock  absorbers  and  bump  stop  energy  dissipation. 
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=*  Initial  conditions  for 
starting  distance  of 
c.g.  from  curb 


Body,  Frame,  Engine 


Contact 

With  Curb  TRWF  TRWR  Decreasing 
2 ' 2 Values 


Negative 

Direction 


Fig.  2 
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3 . 3 Impact  Force 


The  wheel/curb  impact  forces  dependent  on  the 
amount  and  rate  of  deformation  caused  by  curb  impact.  The  impact 
force  is  modeled  using  two  linear  elastic  regions  and  one  region 
of  plastic  deformation  (Fig.  2) . If  the  vehicle  rebounds,  the 
work  done  by  the  deformation  force  across  the  first  and  third 
regions  is  recoverable,  but  the  work  done  across  the  second 
region  is  lost.  The  impact  forces  are  calculated  separately  for 
the  front  and  rear  wheel/suspension  systems. 


3 . 4 Ecfuations  of  Motion 


The  sprung  and  unsprung  masses  are  analyzed  separately  for 
six  of  the  model's  eight  degrees  of  freedom.  The  equations  of 
motion  for  each  masses'  roll  angle  ( <l>s ' lateral  position 
(Ys'  Yu^  vertical  position  (Zg,  Z^)  were  obtained  by  applying 
Newton's  second  law  of  motion  to  the  sprung  mass  (Fig.  3)  and 
unsprung  mass  (Fig.  4) . 

The  forces  which  act  on  the  sprung  mass  include  gravity, 
the  two  forces  of  the  suspension  springs  and  dampers  (Fg^  and 
Fsr) , and  a pin  force  at  the  bushing  which  connects  the  sprung 
and  unsprung  masses  (point  F^  in  Fig.  3 and  4).  The  sprung  mass 
roll  moment  equation  follows: 
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h = ( f"SL  - Fsr)  • s + Fi  • hp^ 


The  equations  of  sprung  mass  lateral  and  vertical  motions  are 
written: 


The  forces  which  act  on  the  unsprung  mass  include  gravity, 
the  suspension  spring  and  damping  forces,  the  bushing  pin  force, 
as  well  as  tire  lateral  friction  forces  (FyL,  ' tire  normal 
reactions  curb  impact  force  (F^)  . The 
equations  of  motion  for  unsprung  mass  roll,  lateral  and  vertical 
motions  are  obtained  respectively: 


= ( ^SL  + ^SR  ) • sin  ( (j)s)  - • cos  ( 


Ms  Zg  - Mg  g = - ( Fsl  + Fsr  ) • cos  ( (l)s)  - F-|  • sin  ( 


z =>  defined  as  zero  when 
at  ground  level 


y 


Fig.  3 
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lu  = ■ ^SL  • s • cos  (A  y + Fsr  • s • cos  (A  4>s)  + F^l  • V - F^r  • x 
- Fyl  * ■ ^YR  * ^1  * *^rc  * ^s^  ' ^ ■*"  *^curb ) 

Mu  Yu  = - ( f^SL  + f^SR ) • sin  ( ^ + • cos  { - FyR  - Fyt  - F^ 

Mu  Zu  - Mu  g = ( Fsl  + Fsr  ) • cos  ( y + F^  • sin  { (!)s)  - F^r  - F^l 

Ai^s  ■ ~^u 

f^ZLF"^  ^ZLR 
*^ZR  = ^ZRF  P'ZRR 

*^CUR8 


Fig.  4 


/ 


The  yaw  equation  of  motion  of  the  vehicle  includes  the 
individual  tire  lateral  forces  (FyLp,  FyLR,  Fyj^p,  ^YRR^  curb 
impact  forces  separately  for  each  axle  (^dF'  ^dR) ' (Fig»  5).  The 
yaw  equation  of  motion  includes  a dynamic  coupling  between  sprung 
mass  roll,  vehicle  pitch,  and  vehicle  yaw,  and  follows: 


I 


= - ( FyLF  + f"YRF  + FdF ) • a + ( FyLR  + FyRR  + F^r  ) 


b + (|)s  0 (^x  ' ^y) 


Vehicle  pitch  is  treated  as  a single  motion  for  the  entire 
vehicle  (Fig.  6) . The  only  forces  affecting  vehicle  pitch  are 
the  individual  normal  tire  reactions  (F^Lp,  ^^ZRF'  ^ZLR'  ^ZRR^  * 
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i 

Fig.  5 

The  pitch  equation  of  motion  includes  a dynamic  coupling  between 
vehicle  yaw,  pitch,  and  sprung  mass  roll  angle: 


ly  0 = ( Fzlf  + Fzrf  ) • a - { Fzlr  + ^ZRR ) * b + r (I^  - Ix) 


(I2-Ix)r<^u  gyroscopic 

moment 

" V 


Fig.  6 
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3 . 5 Modifications  to  STI ' s Tripped  Rollover  Model 


An  inaccuracy  was  discovered  in  an  assumption  made  by  STI 
in  their  tripped  rollover  model.  Several  small  angle  assump- 
tions were  made  affecting  the  sliding  distance  of  the  vehicle 
prior  to  curb  impact  thereby  affecting  the  curb  impact  velocity. 
For  example,  the  distance  that  the  vehicle  skids  before  impact 
should  be  equal  to  the  distance  between  the  vehicle  and  the  curb 
measured  in  the  direction  of  vehicle  motion.  Under  STI ' s 
reference  system,  the  distance  the  vehicle  skids  before  impact 
is  always  measured  perpendicular  to  the  curb,  which  is  a shorter 
distance  in  cases  with  a yaw  angle.  While  STI ' s model  was 
reasonably  accurate  at  small  heading  angles,  it  was  deemed 
necessary  to  add  an  non-inertial  reference  system  (NIRS)  to  the 
model  (Fig.  7) . The  addition  of  the  NIRS  ensures  that  the 
distance  the  vehicle  skids  before  impact  is  perpendicular  to  the 
vehicle  (i.e.  in  the  direction  of  lateral  velocity). 

For  the  purposes  of  sensitivity  analysis  the  STI  simulation 
was  rearranged  so  the  absolute  position  of  the  vehicle  relative 
to  the  curb  was  measured  from  the  front  right  tire.  This  was 
done  to  ensure  that  the  time  of  impact  would  be  independent  of 
vehicle  track  width.  If  the  absolute  position  of  the  vehicle  was 
measured  from  the  vehicle  center  of  gravity  position  (as  it  was 
in  the  original  simulation)  then  the  distance  traveled  prior  to 
impact  by  a vehicle  with  a narrow  track  width  would  be  longer 
than  that  traveled  by  a vehicle  with  a wide  track  width.  Conse- 
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quently  the  two  vehicles  would  hit  the  curb  at  slightly  different 
times,  skewing  the  sensitivity  results.  Measuring  the  vehicles 
absolute  position  from  the  front  right  tire  guarantees  that  the 
distance  traveled  prior  to  impact  is  independent  of  the  vehicle 
track  width,  which  in  turn  guarantees  identical  impact  times. 


3 . 6 Validation  of  Fortran  Subroutine 

In  order  to  perform  a sensitivity  analysis  on  STI ' s tripped 
rollover  model  it  was  necessary  to  translate  STI ' s Basic  program 
into  an  equivalent  Fortran  subroutine.  The  Fortran  version  of 
STI ' s model  was  verified  by  comparing  vehicle  response  and  energy 
using  several  sets  of  initial  conditions.  The  vehicle  data  used 
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throughout  this  project  represents  a 1975  Chevy  Nova  and  is 
listed  in  Appendix  A.  The  system  response  obtained  from  the 
original  STI  simulation  and  the  Fortran  subroutine  were  nearly 
identical  over  the  relevant  time  period  (Fig.  8 - 9) . Small 
differences  appeared  towards  the  end  of  the  rollover  simulation 
when  vehicle  rollover  was  nearly  complete.  The  energy  analysis 
also  had  comparable  results  (Fig.  10  - 11) . We  attribute  any 
discrepancies  between  the  results  obtained  using  the  STI ' s Basic 
simulation  and  UMC's  Fortran  version  to  differences  in  the  met- 
hods used  to  integrate  the  differential  equations  of  motion.  The 
model  was  validated  using  STI ' s reference  system  which  includes 
the  inaccuracies  mentioned  in  section  3.5,  however,  the  change  in 
the  reference  system  should  not  affect  the  system  response  when 
the  heading  angle  was  zero.  Graphs  showing  system  response 
obtained  using  the  Fortran  model  and  1975  Chevy  Nova  data  are 
shown  in  Appendix  C.  The  initial  conditions  (listed  in  Appendix 
B)  represent  a vehicle  which  is  sliding  sideways  toward  a curb  at 
a speed  of  23.075  Ft/S.  The  front  wheel  is  initially  located  7.5 
Ft.  away  from  the  curb. 


SPRUNG  MASS  ROLL  VELOCITY 


(s/peu)  AiioonaA  moH 


16 


YU  D0T=23.075  Ft/S  YO=7 . 5 Ft  H=.0001 


UNSPRUNG  MASS  ROLL  VELOCITY 


(s/PBu)  A1I00R3A  nnoti 


YU  D0T=23.075  Ft/S  Y0=7 . 5 Ft  H=.0001 


RP  ENERGY  RESERVE 


(B^S/5^'4d  aniS)  A0H3N3 
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Ft . ; YU  DOT  = 23.075  Ft . /S 


in 

r- 

II 

o 

> 


0ms) 


A0b!3N3 


19 


Ft . : YU  DOT  = 23.075  Ft . /S 


4 . 0 SENSITIVITY  METHODS  EMPLOYED  IN  ANALYSIS 


A vehicle's  dynamic  system  and  its  subsystems  (suspensions, 
brakes,  steering,  ect.)  are  comprised  of  many  elements  whose 
characteristics  directly  or  indirectly  affect  its  behavior  in 
critical  accident  situations  such  as  those  encountered  during 
vehicle  rollover.  An  understanding  of  how  these  characteristics 
affect  a vehicle's  dynamic  behavior  is  an  essential  part  of 
design  process.  Unfortunately,  the  influence  of  system  para- 
meters on  a vehicle's  response  in  a particular  maneuver  is  not 
always  immediately  evident.  However,  it  is  possible  to  efficien- 
tly determine  the  influence  of  all  parameters  on  system  response 
using  sensitivity  methods,  making  them  very  useful  in  the 
analysis  and  synthesis  of  complex  mechanical  systems. 

Sensitivity  analysis  can  be  carried  out  in  both  the  time  and 
frequency  domains.  When  the  dynamic  characteristics  of  a vehicle 
system  are  given  in  the  form  of  amplitude/phase  frequency  charac- 
teristics or  equivalent  frequency  transfer  functions  then  sensi- 
tivity functions  in  the  frequency  domain  are  the  most  useful.  If 
a vehicle/driver/environment  system  is  given  by  differential 
equations  of  motion  and  the  transient  response  is  of  interest 
then  sensitivity  analysis  in  the  time  domain  is  the  most 
appropriate . 

As  part  of  the  research  sponsored  by  NHTSA  the  University  of 
Missouri-Columbia  has  developed  several  general  purpose  algor- 
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ithms  which  can  be  used  to  calculate  the  sensitivity  functions  of 
mechanical  systems.  Complete  descriptions  of  these  sensitivity 
algorithms  and  the  methods  developed  by  UMC  are  given  in 
references  [7]  and  [8].  The  algorithms  are  capable  of  analyzing 
virtually  any  system,  regardless  of  its  complexity,  and  permit 
analysis  of  nonlinear  transient  mechanical  models  in  the  time 
domain,  nonlinear  steady  state  vehicle  models,  as  well  as  linear 
mechanical  models  in  the  frequency  domain.  The  algorithms  allow 
investigation  of  six  different  types  of  sensitivity  functions 
including  first  order  standard,  first  order  percentage,  first 
order  logarithmic,  second  order  standard,  the  logarithmic  general 
sensitivity  measure,  and  the  percentage  general  sensitivity 
measure.  Deciding  which  type  of  sensitivity  functions  should  be 
used  depends  on  the  objectives  of  the  user.  A brief  description 
of  each  type  of  sensitivity  function  is  given  below. 

4 . 1 First  Order  Standard  Sensitivity  Function 

The  first  order  standard  sensitivity  function  represents  the 
partial  derivative  of  a system  variable  taken  with  respect  to  a 
particular  system  parameter.  This  type  of  sensitivity  function 
is  useful  for  determining  the  influence  of  a particular  parameter 
on  a particular  system  variable.  For  example,  if  a designer 
wished  to  investigate  the  influence  of  wheelbase  on  vehicle  yaw 
velocity  then  first  order  standard  sensitivity  functions  could  be 
used. 
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4 . 2 First  Order  Percentage  Sensitivity  Function 


The  first  order  percentage  sensitivity  function  represents 
the  variable  change  caused  by  a parameter  change  of  some  user 
specified  percentage.  Percentage  sensitivity  functions  are 
useful  for  comparing  the  influence  of  various  parameters  on  a 
particular  system  variable.  For  example,  if  a designer  wanted  to 
know  if  a 2 percent  change  in  wheelbase  effected  yaw  angular 
velocity  more  than  a 2 percent  change  in  track  width  the 
percentage  sensitivity  function  would  show  the  change  in  yaw 
velocity  associated  with  each  parameter  change.  Because 
percentage  sensitivity  functions  have  units  which  are  identical 
to  that  of  the  system  variable,  they  can  be  used  to  compare  the 
effects  of  various  parameters  on  a particular  variable. 

4 . 3 First  Order  Logarithmic  Sensitivity  Function 

The  first  order  logarithmic  sensitivity  function  is  the  only 
dimensionless  sensitivity  function  available,  making  it  useful  to 
compare  and  rank  the  influence  of  various  parameters  on  system 
variables.  Users  should  be  aware  that  if  the  variable  value 
approaches  zero  then  the  logarithmic  sensitivity  functions  can 
become  very  large.  It  should  also  be  noted  that  logarithmic 
sensitivity  is  dependent  on  both  system  sensitivity  and  system 
response,  which  can  sometimes  make  interpretation  of  the  time 
behavior  of  a logarithmic  sensitivity  function  difficult.  How- 
ever, logarithmic  sensitivity  functions  are  useful  in  comparing 
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the  order  of  influence  of  various  system  parameters  at  a particu- 
lar time. 

4 . 4 Second  Order  Standard  Sensitivity  Function 

Second  order  sensitivity  functions  represent  the  second 
partial  derivative  of  a system  variable  taken  with  respect  to  a 
system  parameter.  The  second  order  sensitivity  function  can  be 
used  to  gain  additional  insight  into  results  obtained  using  first 
order  standard  sensitivity  functions.  When  the  system  exhibits 
similar  sensitivity  to  a number  of  parameters  then  second  order 
sensitivity  functions  are  particularly  useful. 

4 . 5 Logarithmic  and  Percentage  General  Sensitivity  Measures 

The  general  sensitivity  measure  determines  the  global  effect 
of  all  selected  parameters  on  a particular  variable.  This  is 
useful  for  analyzing  the  combined  effects  of  related  parameters. 
For  example,  one  could  compare  the  combined  effects  of  parameters 
describing  the  system's  stiffness  with  parameters  describing  the 
system's  mass  properties.  The  general  sensitivity  measure  can  be 
based  on  either  percentage  sensitivity  or  logarithmic  sensitivity 
functions  since  these  types  possess  compatible  dimensions. 

The  different  types  of  sensitivity  functions  can  be  used 
when  carrying  out  sensitivity  analysis.  For  certain  situations 
the  most  recommended  sensitivity  functions  are  given  below. 
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Type  of  Sensitivity  Function 


Application 


First  order  standard 


one  parameter  on  one  variable 


Percentage  and  logarithmic 
First  order  logarithmic 
General  sens,  measure  (perc.) 

General  sens,  measure  (log.) 
Second  order  standard 


- many  parameters  on  one  variable 

- one  parameter  on  several  variables 

~ various  groups  of  parameters 
(geom. , mass,  stiff.,  ect.)  on 
one  variable 

- global  effect  of  many  parameters 
(in  groups)  on  several  variables 

~ when  system  exhibits  similar 
sensitivity  to  a number  of 
parameters 


4 . 6 Numerical  Methods  Used  in  Sensitivity  Analysis 

Because  the  main  objective  of  this  project  was  to  perform  a 
sensitivity  analysis  on  STI's  tripped  rollover  simulation  and 
compare  the  importance  of  system  parameters  on  rollover  behavior 
it  was  decided  that  percentage  sensitivity  functions  best  suited 
the  purposes  of  this  project.  Logarithmic  sensitivity  was  deemed 
inappropriate  since  its  dependence  on  system  response  could  comp- 
licate the  analysis  of  sensitivity  results. 

Time  domain  sensitivity  functions  are  computed  by  integra- 
ting the  mechanical  system's  differential  equations  of  motion 
using  different  sets  of  parameter  values.  The  equations  of 
motion  are  first  integrated  using  nominal  parameter  values.  The 
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parameter  of  interest  is  then  increased  by  some  user  specified 
percentage  (usually  1%)  and  the  equations  of  motion  are 
integrated  again  to  obtain  a varied  system  response.  The  nominal 
parameter  value  is  then  reduced  by  the  user  specified  percentage 
and  the  equations  of  motion  are  integrated  again  to  obtain  a 
second  varied  system  response.  Once  the  system  response  has  been 
determined  using  the  different  sets  of  parameter  values  then  the 
first  order  standard  sensitivity  and  second  order  standard 
sensitivity  functions  can  be  computed  by  taking  the  first  and 
second  partial  derivatives  of  the  variable  of  interest  with 
respect  to  the  parameter  using  a three  point  central  difference 
method.  The  sensitivity  algorithms  are  also  capable  of  using 
either  5 or  7 point  central  difference  method  to  compute  the 
partial  derivatives,  however  the  sensitivity  algorithm  must 
carry  out  additional  model  integrations  to  use  these  methods. 
First  order  logarithmic  sensitivity  is  computed  by  multiplying 
the  first  order  standard  sensitivity  function  by  the  nominal 
parameter  value  and  then  dividing  it  by  the  variable  value 
obtained  using  the  nominal  parameter  set.  Percentage  sensitivity 
is  computed  by  finding  the  difference  in  the  variable  caused  by 
the  user  specified  parameter  change.  Table  1 shows  the 
differences  between  the  various  sensitivity  methods: 
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Table  1. 


Sensitivity 

Type 

Analytical 

Expression 

Numerical 

Expression 

Unit 

1st  Order  Standard 

dv 

0P 

Av 

Ap 

V 

P 

1st  Order  Percentage 

6v  , 

Av 

V 

1st  Order  Logarithmic 

dV  P 

ap  V 

Av  p 
Ap  V 

None 

2nd  Order  Standard 

d^v 

0 p2 

Ap 

V 

p2 

General  Sensitivity 
Measure 

■ n 

. i=l 

-L 

2 (Where  Sj^^  is 
sensitivity 
shown  above) 

any  of  the 
functions 

In  the  case  of  secondary  parameters  (parameters  whose  values 
do  not  remain  constant  e.g.  tire  forces)  the  central  difference 
expressions  cannot  be  used  since  they  assume  that  the  distance 
between  varied  parameter  values  is  constant.  Therefore  under 
these  conditions  the  sensitivity  algorithms  use  a polynomial  fit 
to  determine  the  first  and  second  order  standard  sensitivity 
functions.  The  program  uses  a 2nd  order  polynomial  when  using 
the  three  point  method,  a 4th  order  polynomial  for  the  five  point 
method  and  a 6th  order  polynomial  for  the  seven  point  method.  A 
more  complete  description  of  the  methods  used  in  the  sensitivity 
algorithms  is  given  in  references  [7]  and  [8]. 
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5.0  KNERGY  ANALYSIS  OF  STI  TRIPPED  ROLLOVER  MODEL 


Vehicle  rollover  is  essentially  a problem  associated  with 
stability.  In  a rollover  situation,  a vehicle  which  generates 
small  roll  angles  will  return  to  a four  wheel  stance  and  is 
considered  stable,  however  at  larger  roll  angles  the  vehicle  can 
tip  over  and  is  considered  unstable.  The  stability  of  a system's 
equilibrium  position  can  be  investigated  based  on  the  change  in 
the  system's  potential  energy  caused  by  a small  perturbation  of 
its  state  variables.  If  this  change  is  positive  then  the  system 
is  stable,  if  negative  then  the  system  is  unstable.  For  example, 
when  a vehicle  has  a roll  angle  of  zero  an  increase  in  the 
vehicle  roll  angle  raises  its  center  of  gravity  position,  thus 
increasing  the  potential  energy  of  the  system,  which  indicates 
that  the  initial  position  of  a vehicle  is  stable.  However,  when 
a vehicle's  roll  angle  equals  the  static  tip  over  angle  then 
increasing  the  roll  angle  would  lower  the  c.g.  position, 
decreasing  the  potential  energy  of  the  system,  making  the  vehicle 
unstable.  While  it  was  deemed  impractical  to  perform  a classical 
stability  analysis  on  STI ' s Tripped  Rollover  Vehicle  Simulation, 
it  was  possible  to  examine  the  kinetic  and  potential  energies  of 
the  vehicle.  This  was  done  in  order  to  gain  insight  into  the 
development  of  a function  which  could  be  used  as  an  indication  of 
vehicle  rollover  propensity.  It  was  believed  that  the  maximum 
allowable  change  in  the  gravitational  potential  energy  from  the 
static  condition  to  the  static  tip  over  position  might  serve  as  a 
threshold  which  the  non-centroidal  rotational  kinetic  energy  must 
surpass  in  order  for  vehicle  rollover  to  occur. 
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For  the  sets  of  initial  conditions  considered  in  this 
report,  vehicle  yawing  and  pitching  angular  velocities  were  small 
and  the  kinetic  energies  associated  with  these  motions  were 
neglected.  Under  this  assumption  the  kinetic  energy  of  general 
plane  motion  (Tgpj^)  represents  the  total  kinetic  energy  of  the 
vehicle  and  is  equal  to  the  sum  of  the  system's  kinetic  energies 
due  to  the  sprung  and  unsprung  masses'  c.g.  velocities  and 
rolling  angular  velocities: 


Tgpm 


1 2 1 
2 
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1 2 
+ 2 ® s 
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— I 
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The  kinetic  energy  of  the  general  plane  motion  was  divided 
into  two  parts.  The  kinetic  energy  of  non-centroidal  rotation 
(Tnr)  is  created  during  vehicle  impact  with  the  curb  and 
represents  the  vehicle's  kinetic  energy  in  its  rotational  motion 
around  the  corner  point  of  the  curb: 

^ 1 T 2 1_  ^ 2 

^NR  “ 2 "^2 


The  kinetic  energy  of  translation  (Tijir^r)  is  the  equal  to  the 
kinetic  energy  of  general  plane  motion  minus  the  kinetic  energy 
of  non-centroidal  rotation: 


^TRAN  = ^QPM  - Tnr 


28 


The  gravitational  potential  energy  of  the  vehicle  (Vq)  was 
calculated  relative  to  its  static  equilibrium  position: 

Vg  = ( Mg  + Mu  ) • g • ( h3  - hQ  Q ) 

where:  ^g^g.  = average  C.G.  Height 

h3  = instantaneous  average  C.G.  Height 

Thus,  for  a vehicle  to  attain  a maximum  value  of  gravitational 
potential  energy  the  vehicle’s  center  of  gravity  must  be  placed 
at  the  highest  possible  point  (the  static  tipover  angle) . 

The  change  in  elastic  potential  energy  stored  in  the  suspen- 
sion springs  and  the  vehicle  tires  was  investigated,  however  it 
was  difficult  to  determine  what  portion  of  this  energy  was  a 
contributing  factor  to  rollover  and  what  portion  was  dissipated. 
In  the  final  analysis,  the  change  in  elastic  potential  energy  was 
not  included. 


5.1  Rollover  Prevention  Energy  Reserve  (RPER) 


It  was  felt  that  in  order  for  a vehicle  to  rollover,  the 
kinetic  energy  of  non-centroidal  rotation  (Tj^r) / created  during 
the  vehicle's  impact  with  the  curb,  must  exceed  the  maximum 
allowable  change  in  gravitational  potential  energy  of  the 
vehicle.  A vehicle's  rollover  prevention  energy  reserve  (RPER) 
was  then  developed  using  these  two  quantities  as  an  indicator  of 
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rollover  stability, 
defined: 


The  rollover  prevention  energy  reserve  is 


RPER  = VcRix  - Tnr 


■ 

1 T 

rmw^ 

2 

( Mg  + Mu  ) » g « 

[A 

^ ^c.g/  + 

2 V 

- hc.G. 

1 2 1_  2 
■ 2 2 


= ( Mg  + M^  ) • g • f V 1 + SRSF^  - ij 

It  2 1_  ^ 2 

“ 2 ^2 

TRW 

where  ; SRSF  = r 

2 • hc.G. 


The  kinetic  energy  of  the  vehicle  is  converted  into  gravita- 
tional potential  energy  when  the  vehicle  rolls.  A vehicle  which 
strikes  the  curb  at  the  critical  rollover  side  velocity  (the 
speed  at  which  the  the  vehicle  barely  rolls  over)  should  have 
almost  no  kinetic  energy  at  the  static  rollover  angle  (the  point 
of  maximum  gravitational  potential  energy) . Also,  due  to  energy 
dissipation,  the  vehicle's  kinetic  energy  of  non-centroidal  rota- 
tion must  exceed  its  maximum  gravitational  potential  energy  by 
the  amount  of  which  is  energy  dissipated  by  tire,  suspension  and 
impact  forces  after  the  instant  when  Tj^j^  reaches  its  maximum 
value  and  before  the  static  rollover  angle  has  been  reached. 
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5.2  Results  of  Energy  Analysis 


The  energy  analysis  of  the  tripped  rollover  model  was 
conducted  at  various  speeds  in  order  to  study  the  results  in 
rollover  and  non-rollover  cases.  The  time  histories  of  various 
energy  components  were  examined  (Fig.  12  - 16)  using  two  sets  of 
initial  conditions  in  which  the  vehicle  did  not  rollover  (initial 
side  velocity  equal  to  17  Ft/S  and  22  Ft/S  ) , as  well  as  three 
sets  which  produced  vehicle  rollover  (23.075  Ft/S,  25  Ft/S,  and 
27  Ft/S) . The  initial  distance  of  the  front  wheel  to  the  curb  was 
equal  to  7.5  Ft.  and  the  initial  heading  angle  of  the  vehicle  was 
equal  to  zero  (resulting  in  simultaneous  impact  of  front  and  rear 
wheels)  in  each  of  the  five  cases.  The  lateral  velocity  of 
23.075  Ft/S  was  designated  as  the  critical  rollover  side 
velocity.  The  results  show  that,  after  curb  impact,  the  kinetic 
energy  of  non-centroidal  rotation  increases  to  values  which  are 
greater  than  the  total  kinetic  energy  of  general  plane  motion, 
resulting  in  values  of  the  translational  kinetic  energy  which  are 
less  than  zero.  The  absence  of  several  dynamic  couplings  in  the 
equations  of  motion  in  STI ' s vehicle  model  is  believed  to  be  the 
cause  of  this  discrepancy. 

The  rollover  prevention  energy  reserve  graphs  for  each  of 
the  five  different  test  speeds  are  shown  in  Fig.  17  - 21.  With  a 
lateral  velocity  of  17  Ft/S,  the  vehicle  came  to  a complete  stop 
prior  to  impact  with  the  curb.  Figure  17  shows  a slight  reduc- 
tion in  the  vehicle's  RPER  between  .8  second  and  1 second  due  to 
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the  change  in  chassis  roll  angle  which  returns  to  its  static 
position  when  the  vehicle  comes  to  rest. 

At  the  lateral  velocity  of  22  Ft/S  (Fig.  18)  the  vehicle 
did  impact  the  curb,  but  did  not  rollover.  The  vehicle  impacted 
the  curb,  rolled  to  a maximum  roll  angle  smaller  than  the 
critical  roll  angle  (the  static  tip  over  angle  of  the  vehicle) , 
and  then  returned  to  a stable  position  on  four  wheels.  The 
chassis  continued  to  roll  after  returning  to  a four  wheel  stance, 
causing  rolling  oscillations  which  were  damped  out  by  the  shock 
absorbers.  Figure  18  shows  a significant  drop  in  the  RPER 
immediately  after  curb  impact,  followed  by  oscillations  which 
were  damped  out  as  the  vehicle  came  to  a resting  position.  In 
this  and  all  other  non-rollover  cases  examined,  the  RPER  remained 
positive  throughout  the  simulation. 

The  first  rollover  case  examined  was  at  the  critical 
rollover  lateral  velocity  (23.075  Ft/S,  Fig.  19).  Immediately 
after  curb  impact,  the  vehicle's  RPER  dropped  sharply  to  a value 
of  approximately  -400  Slug  Ft^/S^ . The  RPER  then  increased  as 
the  vehicle's  roll  angle  approached  the  static  tip  over  value. 
At  this  point  the  roll  angular  velocity  of  the  vehicle  was  nearly 
zero.  The  RPER  then  drops  again  once  the  vehicle  passes  through 
the  static  tip  over  position  and  rolls  over. 

The  next  two  cases  examine  the  RPER  using  initial  lateral 
velocities  of  25  Ft/S  and  27  Ft/S  (Fig.  20  and  21) . The  figures 
show  that  the  higher  impact  speeds  result  in  larger  drops  in 
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vehicle  RPER. 


At  an  initial  speed  of  25  Ft/S  the  vehicle's  RPER 
dropped  to  a value  of  -10000  Slug  Ft^/  S^  immediately  after 
impact.  The  vehicle's  RPER  then  increases  to  nearly  zero,  but 
then  falls  again  as  the  vehicle  rolls  over.  At  an  initial  speed 
of  27  Ft/S  the  decrease  in  vehicle  RPER  is  much  greater  (Fig. 
21).  Immediately  after  impact  it's  RPER  drops  to  approximately 
-25000  Slug  Ft^/  S^  and  then  increases  to  -16000  Slug  Ft^  / S^ 
before  falling  again  as  the  vehicle  rolls  over. 

Figure  22  shows  vehicle  RPER  at  each  of  the  five  initial 
speeds.  The  figure  clearly  illustrates  that  the  minimum  value  of 
RPER  is  positive  at  speeds  which  are  smaller  than  the  critical 
rollover  velocity  and  negative  at  speeds  greater  than  the  crit- 
ical rollover  velocity,  which  confirms  that  RPER  is  a very  useful 
indicator  of  vehicle  rollover  propensity  in  tripped  rollover 
maneuvers  with  simultaneous  front  and  rear  wheel  impact. 
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ENERGY  ANALYSIS  OF  STI  MODEL 


(5^S/5„q-d  OnnS)  A9H3N3 
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Ft . : YU  DOT  = 17  Ft/S 


ENERGY  ANALYSIS  OF  STI  MODEL 


(2^s/s^^d  onns)  aousns 
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Ft.:  YU  DOT  = 22  Ft/S 


ENERGY  ANALYSIS  OF  STI  MODEL 


onis)  A9H3N3 
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Ft.:  YU  DOT  = 23.075  Ft/S 
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Ft . ; YU  DOT  = 27  Ft/S 


ROLLOVER  PREVENTION  ENERGY  RESERVE 
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ROLLOVER  PREVENTION  ENERGY  RESERVE 


(E^S/E„Td  OniS)  A0H3N3 
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ROLLOVER  PREVENTION  ENERGY  RESERVE 
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ROLLOVER  PREVENTION  ENERGY  RESERVE 


o 


>- 


00 

T 1 

1 

1 

I 

1 

1 

T 1 

I 

1 

1 

• 

a 

o 

o 

o 

o 

O 

o 

G 

o 

G 

o 

o 

o 

o 

o 

o 

G 

□ 

o 

o 

o 

o 

o 

o 

G 

G 

o 

o 

o 

o 

o 

C 

G 

o 

o 

o 

o 

o 

C 

in 

in 

o 

in 

o 

in 

c 

1 

cu 

n> 

10 

1 

i 

1 

1 

1 

Onns)  A0H3N3 


42 


YU  DOT  = 25  Ft/S 


ROLLOVER  PREVENTION  ENERGY  RESERVE 
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ROLLOVER  PREVENTION  ENERGY  RESERVE 
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6.0  SENSITIVITY  RESULTS 


The  percentage  sensitivity  functions  of  RPER  were  investi- 
gated for  each  primary  (constant)  vehicle  parameter.  All  35 
vehicle  parameters  were  split  into  6 groups  in  order  to  facili- 
tate the  interpretation  of  results.  The  first  two  parameter 
groups  represent  geometrical  parameters  such  as  vehicle  track 
width,  wheelbase,  distance  from  the  axle  to  unsprung  mass  c.g. 
position,  ect.  The  third  group  represents  the  vehicle  mass  and 
inertia  parameters,  the  fourth  stiffness  parameters,  the  fifth 
damping,  and  the  sixth  group  consists  of  parameters  associated 
with  forces  such  as  the  coefficient  of  friction,  force  of  plastic 
suspension  deformation,  ect.  A list  and  description  of  the 
parameters  contained  in  each  group  is  shown  in  Appendix  E. 

The  percentage  sensitivity  functions  of  RPER  were  computed 
for  primary  parameters  using  initial  lateral  velocities  of  22, 
23.075  and  25  Ft/S  and  are  shown  in  Figures  23-40.  In  all  of 
these  cases  the  initial  distance  from  the  front  wheel  to  the  curb 
was  equal  to  7.5  Ft.  and  the  initial  heading  angle  was  zero. 
Vehicle  impact  with  the  curb  occurs  0.3  and  0.4  seconds  depending 
on  the  initial  lateral  velocity.  It  should  be  noted  that  the 
results  obtained  using  the  lateral  speed  of  25  Ft/S  are  plotted 
on  a scale  which  is  different  than  that  used  for  the  22  and 
23.075  Ft/S  cases. 
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6 . 1 Geometric  Parameter  Set  #1 


Figures  23  ~ 25  show  the  influence  of  each  parameter  in 
geometric  set  #1  on  vehicle  RPER  at  each  of  the  three  speeds.  At 
an  initial  side  velocity  of  22  Ft/S  Figure  23  shows  that  the 
track  width  (TRW)  and  the  distance  from  the  axle  to  the  static 
sprung  mass  c.g.  position  (hg)  are  the  most  influential  elements 
of  geometrical  parameter  set  #1  and  exhibit  opposite  influence  on 
vehicle  RPER.  A positive  value  of  sensitivity  taken  with  respect 
to  TRW  indicates  that  increasing  a vehicle's  track  width  in- 
creases its  RPER,  while  a negative  value  of  sensitivity  taken 
with  respect  to  hg  indicates  that  increasing  the  height  of  the 
sprung  mass  c.g.  position  above  the  axle  centerline  decreases 
vehicle  RPER.  The  line  labeled  SRSF  represents  the  static 
rollover  stability  factor  which  is  equal  to  the  half  of  the  track 
width  divided  by  the  center  of  gravity  height.  SRSF  is  a secon- 
dary parameter  which  was  substituted  into  the  model  where  the 
ratio  of  track  width  to  c.g.  height  was  found.  The  sensitivity 
of  RPER  to  this  secondary  parameter  was  found  to  be  a constant 
positive  value  which  was  independent  on  the  initial  velocity. 
These  results  suggest  that  a vehicle  RPER  can  be  increased  by 
increasing  its  static  rollover  stability  factor. 

Since  the  initial  side  velocity  of  22  Ft/S  did  not  result  in 
rollover  the  vehicle  rocks  back  and  forth  after  impact.  Because 
of  this  periodic  motion,  the  sensitivity  of  parameters  TRW  and  hg 
oscillate  as  well,  changing  signs  several  times  during  the  last 
half  of  the  sensitivity  run. 


46 


Overall,  Figure  23  shows  that  a vehicle's  RPER  can  be  in- 
creased by  either  increasing  its  track  width,  or  by  reducing  the 
static  height  of  the  sprung  mass  c.g.  above  the  axle  center  line, 
which  makes  sense  from  a practical  viewpoint.  The  figure  shows 
that  parameters  h^,  hj^,  L,  and  a do  not  appear  to  be  influential 
in  this  case. 

Figures  24  and  25  show  the  percentage  sensitivity  functions 
for  vehicle/curb  impacts  with  initial  lateral  velocities  of 
23.075  Ft/S  and  25  Ft/S,  speeds  which  result  in  vehicle  rollover. 
These  figures  illustrate  that  track  width  (TRW)  and  static  height 
of  the  sprung  mass  above  the  axle  centerline  (hg)  are  still  the 
most  influential  parameters  among  those  shown  in  the  figures,  and 
that  vehicle  RPER  can  be  increased  by  increasing  TRW  or  reducing 

^s- 


Because  the  initial  side  velocities  of  23.075  and  25  Ft/S 
result  in  roll ever  the  vehicle  does  not  rock  back  and  forth  after 
impact  as  it  did  using  an  initial  velocity  of  22  Ft/S.  For  this 
reason  the  percentage  sensitivity  functions  do  not  oscillate  in 
the  rollover  cases. 


6 . 2 Geometric  Parameter  Set  #2 

Figures  26  - 28  illustrate  the  influence  of  geometric  para- 
meter set  #2  on  vehicle  RPER  using  percentage  sensitivity  fun- 
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ctions  at  the  different  initial  lateral  velocities.  Each  of  the 
figures  show  that  the  initial  distance  from  the  curb  is  the 
most  influential  parameter.  This  is  because  changing  the  initial 
distance  from  the  curb  changes  the  time  of  impact  with  the  curb 
as  well  as  speed  of  impact.  This  results  in  the  spikes  seen  in 
Figures  26  - 28.  Since  the  initial  distance  from  the  curb  is  not 
a vehicle  design  parameter  but  rather  a parameter  dependent  on 
circumstance  its  importance  is  not  relevant  and  can  be  ignored. 
With  this  in  mind.  Figure  26  (obtained  using  a speed  of  22  Ft/S) 
indicates  that  the  tire  radius  Tj,  is  very  important.  This 
finding  is  significant  because  an  increase  in  tire  radius 
increases  the  c.g.  height  of  the  entire  vehicle.  The  figure 
shows  that  increasing  the  tire  radius  decreases  a vehicle's  RPER. 
The  figure  also  shows  that  the  undeflected  spring  length  SPRLNG 
is  important  and  that  increasing  a vehicle's  undeflected  spring 
length  decreases  its  RPER.  This  makes  sense  because  in  STI ' s 
model  an  increase  in  SPRLNG  raises  the  c.g.  height  of  the  sprung 
mass.  As  before,  the  oscillation  in  the  sensitivity  functions 
which  appear  after  1 second  are  due  to  the  periodic  roll  motion 
which  occurs  in  the  non-rollover  cases.  Figure  26  also  suggests 
that  vehicle  RPER  can  be  raised  by  increasing  the  height  of  the 
curb. 


The  influence  of  elements  making  up  geometric  parameter  set 
#2  on  vehicle  RPER  at  the  rollover  speeds  of  23.075  and  25  Ft/S 
are  shown  in  Figures  27  and  28.  These  figures  also  reveal  that  a 
vehicle's  RPER  can  be  increased  by  increasing  the  height  of  the 
curb,  decreasing  the  tire  radius,  or  decreasing  the  undeflected 
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spring  length.  The  parameters  s and  ABSDEF  appeared  to  be 
relatively  unimportant. 


6 . 3 Mass  Parameter  Set 

Figures  29  - 31  show  the  influence  of  parameters  associated 
with  vehicle  mass  properties  on  vehicle  RPER.  When  these  figures 
are  compared  with  Figures  25  - 28  it  is  evident  that,  as  a whole, 
the  mass  parameters  exhibit  less  influence  on  the  system  response 
than  geometric  parameters  since  the  magnitudes  of  the  sensitivity 
functions  for  the  mass  parameters  are  much  smaller  up  until 
vehicle  rollover  occurs. 

Figure  29  shows  the  percentage  sensitivity  functions  of 
vehicle  RPER  taken  with  respect  to  each  mass  parameter  using  an 
initial  lateral  velocity  of  22  Ft/S.  The  figure  shows  that 
vehicle  RPER  can  be  raised  by  increasing  vehicle  sprung  mass  and, 
to  a lesser  extent,  vehicle  unsprung  mass.  The  periodic  rolling 
motion  common  to  the  non-rollover  case  causes  an  oscillation  in 
the  sprung  mass  sensitivity  of  vehicle  RPER  at  1.2  seconds. 
Parameters  associated  with  vehicle  inertia  exhibit  virtually  no 
sensitivity  until  the  point  of  impact.  This  is  reasonable  since 
the  angular  velocity  of  the  vehicle  is  nearly  zero  up  until 
impact  with  the  curb.  After  impact  the  sensitivity  values  taken 
with  respect  to  inertial  parameters  remain  smaller  than  those  for 
vehicle  mass. 
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The  sensitivity  results  obtained  at  the  critical  rollover 
speed  of  23.075  Ft/S  (Fig.  30)  reveal  that,  initially,  vehicle 
RPER  can  be  increased  by  increasing  sprung  and  unsprung  masses. 
Later  however,  the  sensitivity  of  the  sprung  mass  changes  sign. 
This  may  be  due  to  the  fact  that  at  small  roll  angles  the  gravi- 
tational force  acting  on  the  sprung  mass  tends  to  return  the 
vehicle  to  a four  wheel  stance  and  is  hence  a stabilizing  force. 
However  at  larger  roll  angles  (angles  at  which  the  line  of  action 
of  the  gravitational  force  on  the  sprung  mass  acts  outside  the 
vehicle’s  track  width)  the  gravitational  force  tends  to  pull  the 
vehicle  over  and  is  a destabilizing  force.  Again  the  sensitivity 
values  taken  with  respect  to  the  inertial  parameters  are  zero  up 
until  the  point  of  impact.  After  impact  the  sensitivity  of  the 
inertial  parameters  remain  small  until  the  sensitivity  functions 
diverge  at  1.3  seconds  at  which  time  they  exhibit  positive 
sensitivity,  indicating  that  vehicle  RPER  can  be  increased  by 
increasing  vehicle  inertia  properties.  This  is  logical  since 
inertial  forces  tend  to  resist  vehicle  motion.  The  divergent 
nature  of  the  curves  is  due  to  large  scale  differences  in  system 
response  caused  by  small  parameter  changes  at  a speed  which  is 
very  near  the  threshold  of  rollover. 

The  sensitivity  of  vehicle  RPER  taken  with  respect  to 
vehicle  mass  parameters  using  an  initial  side  velocity  of  25  Ft/S 
(Fig.  31)  show  trends  which  are  similar  to  those  obtained  using 
initial  side  velocities  of  22  and  23.075  Ft/S.  Again,  vehicle 
RPER  can  be  increased,  at  least  initially,  by  increasing  sprung 
and  unsprung  masses.  After  impact  the  sensitivity  of  the  sprung 
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mass  parameter  changes  sign.  The  sensitivity  of  vehicle  RPER  to 
changes  in  inertial  parameters  remain  zero  up  until  the  time  of 
impact,  after  which  time  the  inertial  parameters  exhibit  small 
positive  sensitivity. 


6.4  Stiffness  Parameter  Set 

Figures  32  - 34  illustrate  the  sensitivity  of  vehicle  RPER 
to  changes  in  vehicle  stiffness  parameters.  A comparison  of 
these  figures  with  those  for  geometrical  and  mass  properties 
(Fig.  23  - 31)  reveals  that  the  variation  of  vehicle  stiffness 
parameters  has  less  of  an  overall  effect  on  vehicle  RPER  than 
variation  of  geometrical  and  mass  parameters. 

Fiaure  32  shows  the  influence  of  stiffness  parameters  on 
vehicle  RPER  using  an  initial  side  velocity  of  22  Ft/S.  The 
figure  shows  that  the  variation  of  stiffness  parameters  have 
little  effect  on  vehicle  RPER.  The  only  two  parameters  which 
exhibit  some  significant  influence  are  suspension  spring  stif- 
fness which  is  negative,  and  the  wheel  compliance  KDl  which 
is  positive. 

At  the  critical  rollover  velocity  of  23.075  Ft/S  the 
only  stiffness  parameter  which  exhibits  a significant  influence 
on  vehicle  RPER  is  suspension  spring  stiffness  K1  (Fig.  33).  The 
divergent  nature  of  the  sensitivity  functions  at  a time  equal  1.3 
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seconds  is  again  due  to  large  scale  differences  in  system 
response  caused  by  small  parameter  changes  at  a speed  which  is 
very  close  to  the  rollover  threshold. 

At  the  initial  lateral  velocity  of  25  Ft/S  the  percentage 
sensitivity  of  vehicle  RPER  to  stiffness  parameter  KD2  reaches  a 
relatively  large  positive  value  just  after  impact,  suggesting 
that  a vehicle's  rollover  propensity  could  be  effected  by  the 
vehicle's  force/deformation  relationship. 


6 . 5 Damping  Parameter  Set 

Figures  35  - 37  illustrate  the  effect  of  damping  parameter 
variation  on  vehicle  RPER.  These  figures  indicate  that,  as  a 
whole,  damping  parameters  have  much  less  influence  on  RPER  than 
geometrical,  mass,  and  stiffness  parameters.  The  shock  absorber 
damping  coefficient  had  the  greatest  amount  of  influence  among 
all  of  the  damping  parameters  examined,  however  its  sensitivity 
value  was  small  and  oscillated  between  positive  and  negative 
values  making  it  difficult  to  draw  any  conclusions  regarding  its 
effect  in  tripped  rollover  situations. 


6 . 6 Force  Parameter  Set 

The  influence  of  force  related  parameters  on  vehicle  RPER  is 
shown  in  Figures  39  - 40.  The  sliding  coefficient  of  friction 
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between  the  tire  and  pavement  is  the  most  influential  parameter 
of  the  set.  This  was  expected  since  this  coefficient  affects  the 
time  at  which  the  vehicle  strikes  the  curb  as  well  as  the  speed 
of  impact.  A higher  coefficient  of  friction  delays  the  time  of 
contact  and  reduces  the  impact  velocity,  resulting  in  a large 
positive  value  of  the  sensitivity  function  at  the  instant  of 
contact. 

At  each  of  the  three  speeds  the  sensitivity  of  deformation 
parameter  POINTl  and  frictional  saturation  speed  K^.  are  quite 
small.  The  deformation  parameter  P0INT2  does  not  exhibit  any 
influence  in  the  case  of  the  initial  side  velocity  of  22  Ft/S, 
suggesting  that  the  impact  deformation  never  reaches  the  third 
deformation  region.  However,  in  the  23.075  and  25  Ft/S  cases 
P0INT2  exhibits  strong  positive  sensitivity  immediately  after 
impact,  suggesting  that  vehicle  RPER  can  be  increased  if  the 
plastic  impact  deformation  is  spread  over  a large  distance.  By 
spreading  the  deformation  over  a large  distance  the  amount  of 
work  done  crushing  the  vehicle  is  increased,  reducing  the  amount 
of  energy  which  can  contribute  to  the  vehicle's  rolling  motion. 

The  constant  force  of  plastic  deformation  WHLDEF  exhibited 
some  influence  in  all  three  cases,  however  it  was  difficult  to 
draw  any  firm  conclusions  regarding  its  effect  on  vehicle 
rollover  propensity. 
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6 . 7 Detailed  Investigation  of  Influential  Parameters 


Once  all  of  the  vehicle's  primary  parameters  had  been 
analyzed,  those  which  exhibited  considerable  influence  on  vehicle 
response  were  analyzed  in  detail.  The  percentage  sensitivity 
functions  of  vehicle  RPER  were  computed  using  5 different  initial 
side  velocities  (17,  22,  23.075,  25,  and  27  Ft/S).  As  before  the 
initial  distance  from  the  front  inside  wheel  to  the  curb  was 
equal  to  7.5  feet  and  the  vehicle's  initial  heading  angle  was 
equal  to  zero.  The  sensitivity  functions  were  computed  for  the 
first  0.8  seconds  of  vehicle  motion  in  order  to  determine  what 
takes  place  during  the  initial  phases  of  the  tripping  maneuver. 
The  six  parameters  which  were  investigated  are,  vehicle  track 
width  (TRW)  , sprung  mass  (Mg)  , tire  radius  (T^^)  , the  horizontal 
distance  from  the  front  axle  to  the  vehicle  c.g.  location  (a) , 
the  static  distance  from  the  sprung  mass  c.g.  position  to  the 
axle  center  line  (hg) , and  the  distance  from  the  axle  center  line 
to  the  lateral  restraining  force  (hj^)  . 

Figure  41  shows  the  sensitivity  of  RPER  taken  with  respect 
to  the  six  parameters  using  an  initial  side  velocity  of  17  Ft/S. 
At  this  speed  the  vehicle  comes  to  a complete  stop  prior  to 
impact  with  the  curb.  The  sensitivity  results  show  that  track 
width  is  the  most  influential  parameter,  followed  by  sprung  mass 
Mg  and  distance  hg.  Parameters  T^,  hj.,  and  a exhibited  very 
little  influence  on  RPER. 
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The  sensitivity  obtained  using  a initial  side  velocity  of  22 
Ft/S  is  shown  on  Figure  42.  The  sensitivity  functions  in  this 
figure  are  identical  to  those  obtained  for  a speed  of  17  Ft/S  up 
until  the  impact  time  of  0.5  seconds.  After  impact  track  width 
becomes  more  important  as  does  hg  and  T^^.  Parameters  a and  h^^ 
exhibit  little  influence  on  the  vehicle  RPER.  The  sensitivity 
results  generated  using  an  initial  velocity  of  23.075  (Fig.  43) 
continue  this  trend. 

At  25  Ft/S  vehicle  track  width  is  still  the  most  influential 
parameter  (Fig.  44) , followed  by  tire  radius  Tj^  and  the  distance 
from  the  static  sprung  mass  c.g.  position  to  the  axle  center  line 
(hg) . Sprung  mass  also  becomes  important,  changing  its  sensiti- 
vity function  from  positive  to  negative  right  after  impact. 
Parameters  a and  h^  still  remain  relatively  unimportant.  These 
trends  are  amplified  in  the  case  using  an  initial  side  velocity 
of  27  Fc/S  (Fig.  45) . 


6 . 8 System  Response  Sensitivity 

The  influence  of  parameters  TRW,  a,  L,  Mg,  1^,  K^,  and  on 
vehicle  system  response  (using  an  initial  side  velocity  of  23.075 
Ft/S)  are  shown  in  Figures  46  - 53.  The  influence  of  these  para- 
meters on  the  sprung  mass  roll  angle  is  shown  in  Figure  46. 
After  curb  impact  (which  occurs  at  approximately  0.45  seconds) 
the  figure  illustrates  that  sprung  mass  roll  angle  can  be  dec- 
reased by  increasing  vehicle  track  width  (TRW) , roll  inertia  of 
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the  sprung  mass  (Ig) , and  the  horizontal  distance  from  the  front 
axle  to  vehicle  c.g.  position  (a) , while  increases  in  vehicle 
sprung  mass  (Mg)  , tire  radius  (Tj_)  and  wheelbase  (L)  will 
increase  the  sprung  mass  roll  angle.  The  suspension  damping 
coefficient  (B^)  , suspension  spring  stiffness  (K^^)  , and  yaw 
moment  of  inertia  (I 2)  each  appear  to  have  little  influence  on 
the  sprung  mass  roll  angle. 

Figures  47  and  48  show  the  influence  of  the  same  parameter 
set  on  sprung  mass  roll  velocity  and  sprung  mass  roll  accelera- 
tion. The  interpretation  of  these  figures  agree  with  that  of 
Figure  46,  however  it  is  made  more  difficult  by  the  discontinuous 
nature  of  the  sensitivity  curves. 

The  influence  of  this  parameter  set  on  the  roll  angle  and 
roll  angular  velocity  of  the  unsprung  mass  are  shown  in  Figures 
49  and  50.  The  conclusions  which  can  be  draw  from  these  figures 
are  quite  similar  to  those  obtained  from  Figure  46.  The  unsprung 
mass  roll  angle  and  angular  velocity  can  be  decreased  by 
increasing  vehicle  track  width  (TRW) , sprung  mass  roll  moment  of 
inertia  (Ig) / and  the  horizontal  distance  from  the  front  axle  to 
the  vehicle  c.g.  position  (a),  whi?e  increases  in  vehicle  sprung 
mass  (Mg)  , tire  radius  (Tj»)  and  wheelbase  (L)  result  in  increases 
in  the  unsprung  mass  roll  angle  and  angular  velocity. 

The  percentage  sensitivity  functions  of  the  sprung  and  un- 
sprung mass  c.g.  heights  to  parameter  changes  are  shown  in  fig- 
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ures  51  & 52.  It  should  be  noted  that  the  reference  system  used 
by  STI's  model  has  its  positive  Z direction  vector  pointing  into 
the  ground.  Therefore  a positive  value  of  percentage  sensitivity 
indicates  that  increasing  the  parameter  value  will  decrease  the 
actual  distance  from  the  center  of  gravity  position  to  the 
ground.  With  this  in  mind  the  figures  show  that  the  distance 
from  the  sprung  mass  c.g.  positions  to  the  ground  can  be  in- 
creased by  increasing  tire  radius  (Tj^)  , spring  stiffness  (K-]^)  , 
and  after  impact  track  width  (TRW) , while  the  distance  from  the 
unsprung  mass  c.g.  position  to  the  ground  can  be  increased  prim- 
arily by  increasing  tire  radius  (T^) . Increasing  the  sprung  mass 
(Mg)  will  decrease  the  distance  from  the  sprung  mass  c.g.  posi- 
tion to  the  ground.  It  appears  that  the  parameters  a,  L,  Ig,  I^, 
and  Bjl  have  little  influence  on  the  c.g.  heights  of  both  sprung 
and  unsprung  masses. 

The  influence  of  the  parameters  on  sprung  mass  vertical 
velocity  is  shown  in  Figure  53.  The  sensitivity  functions  taken 
with  respect  to  parameters  T^^  and  Mg  oscillate  between  positive 
and  negative  values  which  makes  the  interpretation  of  results 
difficult.  Parameters  a,  L,  Ig,  1^,  and  all  have  little 
influence  on  the  vertical  velocity  of  the  sprung  mass. 


6 . 9 Sensitivity  of  RPER  to  System  Response 

The  final  sensitivity  functions  presented  in  this  report 
investigate  the  influence  of  system  state  variables  on  vehicle 
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RPER.  Figures  54  - 58  show  how  the  variation  of  Zg,  Zg,  Z^, 
Yg,  and  affects  RPER  at  the  initial  lateral  velocities  of  17, 
22,  23.075,  25,  and  27  Ft/S.  At  each  speed  the  figures  indicate 
that  the  state  variables  Zg,  Z^,  Yg,  and  Y^  have  little  influence 
on  vehicle  RPER,  however,  the  figures  do  show  that  vehicle  RPER 
can  be  increased  by  reducing  the  sprung  and  unsprung  center  of 
gravity  heights  Zg  and  Z^. 
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7 . 0  CONCLUSIONS  AND  RECOMMENDATIONS 


Based  on  the  results  of  the  sensitivity  analysis  the  follow- 
ing conclusions  were  reached. 

7 . 1  Rollover  Prevention  Enercry  Reserve 

The  stability  based  concept  of  rollover  prevention  energy 
reserve  (RPER)  is  a useful  tool  for  examining  the  influence  of 
vehicle  parameters  in  tripped  rollover  maneuvers.  The  results 
reveal  that  the  minimum  value  of  RPER  approaches  zero  at  the 
threshold  of  rollover,  and  if  the  minimum  falls  below  zero  then 
vehicle  rollover  will  occur. 


7 . 2  Sensitivity  Results 

In  general,  the  results  obtained  using  the  percentage  sens- 
itivity functions  are  reasonable  and  agree  with  the  existing 
information  and  opinions  regarding  rollover  behavior  and 
propensity.  While  this  in  itself  does  not  validate  the 
sensitivity  methods  used  in  this  report  it  certainly  lends  them 
credence. 

7.2.1  Geometric  Parameters 

The  sensitivity  results  show  that  vehicle  geometric 
parameters  are  by  far  the  most  influential  characteristics  in  a 
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tripped  rollover  situation.  The  results  indicate  that,  in  trip- 
ped rollover  maneuvers,  vehicle  stability  can  be  best  enhanced  by 
increasing  vehicle  track  width  or  reducing  vehicle  c.g.  height. 
It  was  found  that  geometric  parameters  commonly  associated  with 
suspension  components,  such  as  the  distance  separating  the 
suspension  springs  (s) , and  the  distance  from  roll  axis  to 
vehicle  axle  (hj.)  , exhibited  very  little  influence  on  vehicle 
RPER. 

7.2.2  Mass  Parameters 

The  results  show  that  parameters  generally  associated  with 
the  vehicle's  mass  properties  are  relatively  important,  though 
not  nearly  as  important  as  parameters  associated  with  vehicle 
geometry. 

The  sensitivity  of  RPER  to  variations  in  mass  parameters 
indicate  that  the  stability  can  be  increased  by  increasing  the 
vehicle  sprung  mass  (Mg) . This  result  may  be  partially  due  to 
the  fact  that  the  a heavier  sprung  mass  compresses  the  suspension 
springs,  thus  lowering  the  c.g.  position  of  vehicle.  The  re- 
sults also  reveal  that  increasing  roll  moments  of  inertia  tends 
to  increase  RPER  by  opposing  vehicle  motion. 

7.2.3  Stiffness  Parameters 

The  influence  of  vehicle  stiffness  parameters  on  vehicle 
RPER  was  found  to  be  quite  small,  indicating  that  suspension 
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spring  stiffness  has  little  influence  on  vehicle  motion  in  trip- 
ped rollover  situations. 

7.2.4  Damping  and  Force  Parameters 

Damping  parameters  such  as  the  damping  coefficients  of  shock 
absorbers  were  found  to  have  very  little  influence  on  vehicle 
stability  in  tripped  rollover  maneuvers. 

Parameters  associated  with  external  forces  acting  on  the 
vehicle  such  as  deformation  forces  caused  by  impact,  and  the 
points  of  discontinuity  in  the  force/deformation  relationship 
were  found  to  have  a significant  influence  on  vehicle  stability. 
By  spreading  the  impact  deformation  over  a large  distance  (thus 
increasing  the  amount  of  vehicle  plastic  deformation)  the  amount 
of  unrecoverable  work  used  to  crush  the  vehicle  is  increased, 
reducing  the  amount  of  energy  which  can  contribute  to  the 
vehicle's  rolling  motion. 


7 . 3 Recommendations 


Accident  statistics  have  shown  that  vehicle  rollover  rep- 
resents a serious  safety  problem  on  the  nations  roads  and  high- 
ways and  that  utility  vehicles  and  small  pickup  trucks  have  the 
highest  probability  of  being  involved  in  a rollover  accident. 
Since  the  sales  of  these  vehicles  have  grown  substantially  in  the 
past  few  years  it  is  likely  that  the  number  of  fatalities  due  to 
rollover  accidents  will  increase  as  well.  Because  of  these  trends 
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there  has  been  a growing  call  for  the  government  to  develop 
legislation  which  will  establish  minimum  stability  standards  for 
passenger  vehicles.  If  legislation  aimed  at  curbing  the  number 
of  rollover  accidents  is  to  be  effective  it  is  imperative  that 
the  influence  of  design  parameters  on  vehicle  rollover  propensity 
be  thoroughly  understood.  With  this  in  mind  we  feel  that  more 
research  into  vehicle  rollover  behavior  is  necessary  and  should 
be  directed  toward  the  following  areas. 

7.3.1  Continued  Investigation  Into  Tripped  Rollover 

This  report  presents  a thorough  investigation  into  the  ef- 
fects of  vehicle  design  parameters  on  tripped  rollover  maneuvers 
involving  simultaneous  impact  of  front  and  rear  wheels.  It 
should  be  noted,  however  that  the  results  presented  in  this  paper 
were  obtained  using  a single  set  of  vehicle  data  and  in  one 
particular  type  of  rollover  accident.  Because  of  the  this  we 
feel  that  the  sensitivity  methods  used  in  this  report  should  be 
applied  to  investigate  several  other  sets  of  vehicle  data  and  in 
different  rollover  situations  using  an  improved  tripped  rollover 
vehicle  model.  The  vehicle  data  sets  should  reflect  the  dif- 
ferent classes  of  vehicles  which  are  currently  in  use.  These 
classes  should  include  small  and  subcompact  passenger  vehicles, 
small  and  standard  pickup  trucks,  and  utility  vehicles. 
Maneuvers  to  be  investigated  should  include  simultaneous  front 
and  rear  wheel  impact  as  well  as  vehicle  oblique  impact  with  the 
curb.  An  improved  vehicle  model  should  be  developed  which  would 
account  for  all  dynamic  couplings  in  the  equations  of  motion  and 
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allow  for  a thorough  accounting  of  energy  transfer  and  its 
dissipation. 

7.3.2  Investigation  Into  Untripped  Rollover  Behavior 

A great  deal  of  the  controversy  surrounding  the  rollover 
issue  involves  the  apparent  propensity  of  certain  utility 
vehicles  to  rollover  in  collision  avoidance  maneuvers  without  the 
help  of  a tripping  mechanism.  The  influence  of  vehicle  design 
parameters  on  rollover  propensity  in  these  types  of  untripped 
maneuvers  may  differ  considerably  from  those  associated  with 
tripped  rollover.  While  this  report  found  that  the  influence  of 
parameters  associated  with  suspension  characteristics  in  tripped 
rollovers  were  small  these  parameters  may  prove  to  have  a 
substantial  effect  on  vehicle  behavior  in  untripped  maneuvers. 
Because  of  this  it  would  be  worthwhile  to  develop  a vehicle  model 
which  would  be  capable  of  predicting  the  occurrence  untripped 
rollover  of  vehicles  in  collision  avoidance  maneuvers.  Once  this 
is  accomplished  a sensitivity  analysis  of  the  model  should  be 
performed  to  determine  the  influence  of  vehicle  design  parameters 
on  rollover  propensity  in  untripped  situations. 
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APPENDIX  A 


Vehicle  Data 

for  1975  Chew  Nova 

a 

4.45 

ft 

L 

9.25 

ft 

TRW 

5.00 

ft 

0.2 

ft 

Kz 

12000.00 

lb/ ft 

80.00 

Ib-sec/ft 

0.75 

— 

Bd 

0.25 

Ib-sec/ft 

KDl 

24000.00 

lb/ ft 

KD2 

96000.00 

lb/ ft 

Mu 

16.7 

lb-s^2/ft 

hr 

1.0 

ft 

Kl 

6900.00 

lb/ ft 

^2 

33000.00 

lb/ ft 

325.00 

Ib-sec/ft 

^2 

33 . 00 

Ib-sec/ft 

330.00 

lb-ft/s^2 

2240.00 

lb-ft/s^2 

iz 

2240.00 

Ib-ft/s  ^2 

1.2 

ft 

lu 

132 . 00 

lb-ft/s^2 

WHLDEF 

4000.00 

lb 

Kv 

1.00 

ft/s 

POINTl 

0.1667 

ft 

A-1 


Vehicle  Data  for  1975  Chew  Nova  (Concluded) 


P0INT2 

0.6667 

ft 

Ms 

105.00 

lb-S^2/ft 

240.00 

lb-ft/s^2 

T 

1.0 

ft 

Mcurb 

0.5 

ft 

s 

1.75 

ft 

K 

200000.00 

lb/ ft 

SPRLNG 

0.6667 

— 

ABSDEF 

0.25 

ft 

B 

1000.00 

Ib-s/ft 

A-2 


APPENDIX  B 


Initial  Conditions  for  Tripped  Rollover  Maneuvers 


Yu 

dot 

varied* 

ft/s 

Zu 

dot 

0.00 

ft/s 

Phi^ 

dot 

0.00 

rad/s 

Ys 

dot 

varied* 

ft/s 

Zs 

dot 

0.00 

ft/s 

Phis 

dot 

1 

0.00 

rad/s 

Theta 

dot 

0.00 

rad/s 

Psi 

dot 

0.00 

rad/s 

Yu 

10.00** 

ft 

Zu 

-1.11 

ft 

Phi^ 

1.1568 

rad 

Ys 

10.00** 

ft 

Zs 

-2.5319 

ft 

Phis 

1.8437 

rad 

Theta 

0.00 

rad 

Psi 

0.00 

rad 

* Initial  lateral  speeds  used  were  17  ft/s,  22  ft/s,  23.075  ft/s, 
25  ft/s,  and  27  ft/s 

**  Vehicle  initial  position  was  7.5  ft  measured  from  the  front 
right  tire,  the  location  of  the  added  non-inertial  reference 
system. 
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APPENDIX  C 


System  Response  of  Vehicle  in  Tripped  Rollover 
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APPENDIX  D;  Fortran  Source  Code  of  STI ' s Tripped 

Vehicle  Rollover  Model 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

c 


c 

c 


c 

c 

c 


***************************************************************** 


* * 

* SUBROUTINE  NAME:  MODEL  * 

* * 

* PURPOSE:  8 DOF  NON-LINEAR  MODEL  BASED  ON  SYSTEMS  * 

* TECHNOLOGY  INC.  ROLL  OVER  MODEL  * 

* * 

* CALLED  BY  SUBROUTINE:  FCN  * 

* * 

* LAST  REVISION  - 4/28/88  BY  - ANDRZEJ  G.  NALECZ  * 

* ALAN  BINDEMANN  * 

* * 

* CALLING  SEQUENCE:  * 

* * 

* CALL  MODEL(N,X, Y, YPRIME,VAR)  * 

* * 

* ARGUMENTS : * 

* * 

* INPUT:  COMMON  BLOCKS  VDAT,  AND  SECP  * 

* * 

* OUTPUT:  INTEGRATION  ARGUMENTS  D(l-16)  * 

* * 


*************************************************************** 


SUBROUTINE  MODEL ( N , X , Y , YPRIME , VAR , FLG 1 ) 
IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 


COMMON/VDAT/ 

A, 

L, 

TRW, 

HU, 

KZ, 

BZ, 

1 

CF, 

BD, 

KDl, 

KD2  , 

MU, 

HR, 

SUSKl, 

SUSK2 , 

2 

Bl, 

B2  , 

IX, 

lY, 

IZ, 

HS, 

lU, 

WHLDEF 

3 

KV, 

POINTl, 

POINT2 , 

MS, 

IS, 

TRAD, 

HCURB , 

s, 

4 

SPRNGK, 

SPRLNG, 

ABSDEF, 

BL, 

YO, 

DUM, 

PCI, 

PC2  , 

5 

PC3  , 

PC4, 

PC5, 

PC6, 

PC7, 

PC8  , 

PC9  , 

PCIO, 

6 PCll 


COMMON/ STATUS/  IP,  INM 

COMMON/SECP/  SRSF,  SPCF,  DAMF,  SPCR,  DAMP,  ZSP,  ZUP,  ZSDP, 
1 ZUDP,  YSDP,  YUDP 

INTEGER*4  N, IP, INM 


1 

2 

3 

4 

5 

6 


REAL* 8 MS,  MU, 

KZ,  KV, 

MAXTRW, 
TWR(50,7) , 
TRWR(50, 7) , 
YRMAX(50, 7) , 
VAR (50) , 


IS,  lU, 

ISA,  L, 

H4(50,7) , 
CRIT(50, 7) , 
YINT2 (50,7) , 
MAXTRWF(50, 7) , 
YPRIME (17) , 


IX,  lY,  IZ, 

KDl,  KD2,  lUA, 

TWF(50,7) , 

TRWF(50,7) , 
YFMAX(50,7) , 
MAXTRWR(50, 7) , 

Y(17) 


(1)  ASSIGN  INITIAL  VALUES 

(2)  UPDATED  PARA  SENT  BACK  TO  MODEL 


D-1 


ooo  ooo  ooo 


c 


CALL  VDATA 


ASSIGN  VARIABLE  NAMES 


YSD 

= 

Y(l) 

ZSD 

= 

Y(2) 

PHISD 

= 

Y(3) 

YUD 

= 

Y(4) 

ZUD 

= 

Y(5) 

PHIUD 

= 

Y(6) 

THED 

= 

Y(7) 

YAWR 

= 

Y(8) 

YS 

= 

Y(9) 

ZS 

= 

Y(10) 

PHIS 

= 

Y(ll) 

YU 

= 

Y(12) 

ZU 

= 

Y(13) 

PHIU 

= 

Y(14) 

THE 

= 

Y(15) 

PS  I 

= 

Y(16) 

YGL 

= 

Y(17) 

CALCULATE  BUMP  STOP  SIZES 


G = 32.17D0 

BSDEF=SPRLNG-MS*G/SUSKl/2 . ODO 
BSSIZL  = BSDEF  - ABSDEF 
BSSIZU  = BSDEF  + ABSDEF 


CALCULATE  INITIAL  VARIABLES  (FIRST  MODEL  EVALUATION  ONLY) 


C 


IF  (X.EQ.O.ODO)  THEN 
CRIT(IP,INM) 
TWF(IP,INM) 
TWR(IP,INM) 
TRWF(IP,INM) 
TRWR(IP,INM) 
MAXTRWF(IP,INM)  = 
MAXTRWR(IP, INM)  = 


DAT  AN  (TRW/ 2 . 0 DO/HU) 

TRW/ 2. ODO 

TRW/ 2. ODO 

TRW/2. ODO 

TRW/ 2. ODO 

TRW/2 . ODO 


^ TRW/ 2. ODO 

PHIU  = (MU*G*CF* (TRAD  + HU) +MS*G*CF* (HS+TRAD) ) /KZ/TRW/TRW 
PHIS  = (MS*G*CF*(HS  - HR) /2 . ODO/SUSKl/S/S)  + PHIU 
ZUP  = ZU  * PC7 
ZS  = ZUP  - (HS  - HU  + BSDEF) 

ZSP  = ZS  * PC6 
YS  = YU  + (ZUP  - ZSP) 

YSDP  = YSD  * PCIO 
YUD  = YSDP 
YUDP  = YUD  * PCll 
ZSDP  = ZSD  * PCS 
ZUDP  = ZUD  * PC9 


* DCOS(PHIS) 
* DTAN(PHIS) 


H4(IP,INM)  = -(MS  * ZSP  + MU  * ZUP)  / (MS  + MU) 
YINT2 (IP, INM)  = WHLDEF  - KD2  * POINT2 
END  IF 


D-2 


ZSP  = zs 

* 

PC6 

ZUP  = zu 

* 

PC7 

ZSDP  = ZSD 

* 

PCS 

ZUDP  = ZUD 

* 

PC9 

YSDP  = YSD 

* 

PCIO 

YUDP  = YUD 

* 

PCll 

c 

C CHECK  FOR  ROLLOVER  OR  POTENTIAL  OVERFLOW  IN  DGEAR 
C 

IF  (INM.Eq.l)  THEN 

IF  (DABS (PHIU) .GT.CRIT(IP, INM) ) THEN 
WRITE (5, 105)  PHIU  * 57.3D0 

105  FORMAT ( IX PHIU  = ' , F7 . 3 , ' DEGREES') 

WRITE (5,*)  'UNSPRUNG  MASS  ROLL  ANGLE  > THAN  CRIT.  ANGLE' 
WRITE (5,*)  'ROLLOVER  ASSUMED' 

STOP 
END  IF 

ELSE  IF  (DABS (PHIU) .GT. 1. 57 08D0)  THEN 
WRITE(5,106)  PHIU  * 57.3D0 

106  FORMAT ( IX ,' PHIU  = ' , F7 . 3 , ' DEGREES  DURING  PARAM.  VARIATION') 
WRITE (5, 107)  IP,INM 

107  FORMAT ( IX, ' IP  = ',12,'  INM=  *,I2) 

WRITE (5,*)  'COMPUTATIONS  STOPPED  TO  PREVENT  OVERFLOW  IN  DGEAR' 
STOP 


END  IF 

SINPSI 

= 

-DSIN(PSI) 

COSPSI 

= 

DCOS (PSI) 

SINPHIU 

= 

DSIN(PHIU) 

COSPHIU 

= 

DCOS (PHIU) 

TANTHE 

= 

DTAN(THE) 

COSTHE 

= 

DCOS (THE) 

COS2 

= 

COSPHIU  * COSPHIU 

COS2THE 

= 

COSTHE  * COSTHE 

B=L-A 

C 

C CALCULATE  CRUSH  DISTANCE  IN  Y DIRECTION 
C 


YF=(YO-YGL)/COSPSI  ! UMC  DEFORMATION 

YR=YF+(A+B) *DTAN(PSI)  ! UMC  DEFORMATION 

C 

C CALCULATE  RATE  OF  DEFORMATION 
C 


C 

C 


YFD=(Y0-YGL) *YAWR*DTAN ( PSI ) /COSPSI-YUDP-YAWR* (A- 
1 TRW/2 . 0D0*DTAN( PSI) ) 

YRD=YFD+ (A+B) *YAWR/COSPSI/COSPSI 


(COSPSI.lt. 0. ODO)  THEN 

! UMC 

DEFORMATION 

YF=-YF 

! UMC 

DEFORMATION 

YR=-YR 

! UMC 

DEFORMATION 

YFD=-YFD 

! UMC 

DEFORMATION 

YRD=-YRD 

! UMC 

DEFORMATION 

1 IF 

! UMC 

DEFORMATION 

LOGIC 

LOGIC 


LOGIC 

LOGIC 

LOGIC 

LOGIC 

LOGIC 

LOGIC 


! UMC  LOGIC 
! UMC  LOGIC 
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LEFT  FRONT  TIRE 

ZT  = (ZUP  - TRW/2. ODO  * SINPHIU) /COSPHIU  + HU  - A * TANTHE 
ZTD  = ZUDP/COSPHIU  - PHIUD  * ( TRW/ 2 . ODO  - ZUP*SINPHIU) /COS2 
1 - A * THED/COS2THE 

ZTLF  = ZT 
ZTDLF=ZTD 

FZT=KZ* (TRAD+ZT) +BZ*ZTD 
IF  (FZT.LE. 0. ODO)  GOTO  467 
GOTO  468 

67  FZT  = O.ODO 

68  CONTINUE 
FZLF  = FZT 
DIST  = HU  + DABS(ZT) 

COSBETA=DCOS (1 . 5707-DABS (PHIU) -DATAN(DIST/ (TRW/ 2 . ODO) ) ) 

ROLL  = DSQRT(DIST*DIST  + TRW*TRW*0 . 25D0)  * COSBETA 
VEL=YUDP+A*YAWR-PHIUD*ROLL 
FYV=VEL/KV 

IF  (VEL.GT.  KV)  FYV  = l.ODO 
IF  (VEL.LT.-KV)  FYV  = -l.ODO 
FYT  = CF  * FZT  * FYV 
FYLF  = FYT 

RIGHT  FRONT  TIRE 

ZT  = ( ZUP  + TRWF(IP,INM) *SINPHIU  ) /COSPHIU  + HU  - A*TANTHE 
ZTD  = ZUDP/COSPHIU  + PHIUD  * (TRWF ( IP , INM)  + ZUP  * SINPHIU)/ 

1 COS2  - A * THED/COS2THE 

ZTRF  = ZT 
ZTDRF  = ZTD 

FZT=KZ* (TRAD+ZT) +BZ*ZTD 
IF  (FZT.LE. 0. ODO)  GOTO  1467 
GOTO  1468 

1467  FZT  = O.ODO 

1468  CONTINUE 
FZRF  = FZT 

DIST  = HU  + DABS(ZT) 

COSBETA  = DCOS(1.5707  - DABS(PHIU)  - DATAN (DIST/TRWF (IP, INM) ) ) 
ROLL  = DSQRT(DIST*DIST  + TRWF (IP, INM) *TRWF (IP, INM) ) * COSBETA 
VEL  = YUDP  + A*YAWR  - PHIUD*ROLL 
FYV=VEL/KV 

IF  (VEL.GT.  KV)  FYV  = l.ODO 
IF  (VEL.LT.-KV)  FYV  = -l.ODO 
FYT  = CF  * FZT  * FYV 
FYRF  = FYT 
C 

C CALCULATE  THE  FRONT  RIGHT  DEFORMATION  FORCE 
C 

C YF=Y0-YU+A*SINPSI-TRW/2.0D0*COSPSI  ! STI  LOGIC 

C YFD=-YUDP+A*YAWR*COSPSI+TRW/2 . 0D0*YAWR*SINPSI  ! STI  LOGIC 

C 

DUMY=YF 

DYMAX=YFMAX ( IP , INM) 

MAXTRW=MAXTRWF ( IP, INM) 
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DTRW=TRWF ( I P , I NM ) 

C 

ABD  = 1.0 

IF  (DUMY.GE. 0. ODO)  GOTO  501 
DTRW  = DUMY+TRW/2 . ODO 
DDUMY  = -DUMY 

C DTRW  = ( (YO-YU+A*SINPSI)/COSPSI) 

C DDUMY  = TRW/2.0D0-DTRW 

IF  (DDUMY. LT.DYMAX)  GOTO  510 
DYMAX  = DDUMY 

IF  ( DDUMY. GT.POINTl)  GOTO  496 
DUMYF  = KDl  * DABS (DDUMY) 

GOTO  518 

496  IF  ( DDUMY. GT.POINT2)  GOTO  499 

DUMYF  = WHLDEF 
GOTO  518 

499  DUMYF  = KD2  * DABS (DDUMY)  + YINT2 ( IP, INM) 

GOTO  518 

501  DUMYF  = O.ODO 

ABD  = O.ODO 
GOTO  521 
C 

C FRONT  CURB  FORCE  AFTER  MAXIMUM  DEFORMATION 
C 

510  IF  (DYMAX. LT.POINT2)  GOTO  513 

DUMYF  = KD2  * DABS (DDUMY)  + YINT2 (IP, INM) 

GOTO  517 

513  IF  (DYMAX. LT.POINTl)  GOTO  516 

DUMYF  = WHLDEF  - KDl  * DABS (DYMAX-DDUMY) 

GOTO  517 

516  DUMYF  = KDl  * DABS (DDUMY) 

517  IF  (DUMYF.lt. 0. ODO)  THEN 

DUMYF  = O.ODO 
ABD  = O.ODO 
END  IF 

518  IF  (DTRW.GT.MAXTRW)  GOTO  520 

MAXTRW  = DTRW 

520  DTRW  = MAXTRW 

521  TRWDUM  = MAXTRW 
C 

YFMAX(IP, INM)  = DYMAX 
TRWF(IP,INM)  = DTRW 
TWF(IP,INM)  = TRWDUM 
MAXTRWF(IP, INM)  = MAXTRW 
SPCF  = DUMYF  * PC2 
DAMF  = BD  * ABD  * YFD  * PC3 
FDF  = SPCF  - DAMF 
C 

C LEFT  REAR  TIRE 
C 

ZT  = (ZUP-TRW/2 . ODO*SINPHIU)/COSPHIU  + HU  + B*TANTHE 
ZTD  = ZUDP/COSPHIU  - PHIUD* (TRW/2 . ODO  - ZUP*SINPHIU)/COS2  + 
1 B * THED/COS2THE 

ZTLR  = ZT 
ZTDLR  = ZTD 


! UMC  DEFORMATION  LOGIC 
! UMC  DEFORMATION  LOGIC 
! STI  DEFORMATION  LOGIC 
! STI  DEFORMATION  LOGIC 
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FZT  = KZ  * (TRAD  + ZT)  + BZ  * ZTD 
IF  (FZT.LE. 0. ODO)  GOTO  2467 
GOTO  2468 

2467  FZT  = O.ODO 

2468  CONTINUE 
FZLR  = FZT 
DIST  = HU  + DABS(ZT) 

COSBETA  = DCOS (1.5707  - DABS(PHIU)  - DATAN (DIST/ (TRW/2 . ODO) ) ) 
ROLL  = DSQRT(DIST*DIST  + TRW*TRW*0 . 25D0)  * COSBETA 
VEL  = YUDP  - B*YAWR  - PHIUD*ROLL 
FYV=VEL/KV 

IF  (VEL.GT.  KV)  FYV  = l.ODO 
IF  (VEL.LT.-KV)  FYV  = -l.ODO 
FYT  = CF  * FZT  * FYV 
FYLR  = FYT 

RIGHT  REAR  TIRE 

ZT  = ( ZUP  + TRWR(IP, INM) *SINPHIU  ) /COSPHIU  + HU  + B*TANTHE 
ZTD  = ZUDP/ COSPHIU  + PHIUD  * (TRWR ( IP , INM)  + ZUP  * SINPHIU)/ 

1 COS2  + B * THED/COS2THE 

ZTRR  = ZT 
ZTDRR  = ZTD 

FZT=KZ* (TRAD+ZT) +BZ*ZTD 
IF  (FZT.LE. 0. ODO)  GOTO  3467 
GOTO  3468 

67  FZT  = O.ODO 

68  CONTINUE 
FZRR  = FZT 
DIST  = HU  + DABS(ZT) 

COSBETA  = DCOS (1.5707  - DABS(PHIU)  - DATAN (DIST/TRWR ( IP, INM) ) ) 
ROLL  = DSQRT(DIST*DIST  + TRWR(IP, INM) *TRWR(IP, INM) ) * COSBETA 
VEL  = YUDP  - B*YAWR  - PHIUD*ROLL 
FYV=VEL/KV 

IF  (VEL.GT.  KV)  FYV  = l.ODO 
IF  (VEL.LT.-KV)  FYV  = -l.ODO 
FYT  = CF  * FZT  * FYV 
FYRR  = FYT 

CALCULATE  THE  RIGHT  REAR  DEFORMATION  FORCE 

YR=Y0-YU-B*SINPSI-TRW/2 . ODO*COSPSI  ! STI  LOGIC 

YRD=-YUDP-B*YAWR*COSPSI+TRW/2.0D0*YAWR*SINPSI  ! STI  LOGIC 

DUMY=YR 

DYMAX=YRMAX ( IP , INM) 

MAXTRW=MAXTRWR ( I P , I NM ) 

DTRW=TRWR ( IP, INM) 

ABD  = 1.0 

IF  (DUMY.GE. O.ODO)  GOTO  701 
DTRW  = DUMY+TRW/2 . ODO 
DDUMY  = -DUMY 

DTRW  = ( (YO-YU-B*SINPSI)/COSPSI) 

DDUMY  = TRW/2 . ODO-DTRW 


! UMC  DEFORMATION  LOGIC 
! UMC  DEFORMATION  LOGIC 
! STI  DEFORMATION  LOGIC 
! STI  DEFORMATION  LOGIC 
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IF  (DDUMY.lt. DYMAX)  GOTO  710 
DYMAX  = DDUMY 

IF  ( DDUMY. GT.POINTl)  GOTO  696 
DUMYF  = KDl  * DABS (DDUMY) 

GOTO  718 

696  IF  ( DDUMY. GT.POINT2)  GOTO  699 

DUMYF  = WHLDEF 
GOTO  718 

699  DUMYF  = KD2  * DABS(DDUMY)  + YINT2 ( IP, INM) 

GOTO  718 

701  DUMYF  = O.ODO 

ABD  = O.ODO 
GOTO  721 
C 

C REAR  CURB  FORCE  AFTER  MAXIMUM  DEFORMATION 
C 

710  IF  (DYMAX. LT.POINT2)  GOTO  713 

DUMYF  = KD2  * DABS (DDUMY)  + YINT2 (IP, INM) 
GOTO  717 

713  IF  (DYMAX.lt. POINT 1)  GOTO  716 

DUMYF  = WHLDEF  - KDl * DABS (DYMAX-DDUMY) 
GOTO  717 

716  DUMYF  = KDl  * DABS (DDUMY) 

717  IF  (DUMYF.lt. 0. ODO)  THEN 

DUMYF  = O.ODO 
ABD  = O.ODO 
END  IF 

718  IF  (DTRW.GT.MAXTRW)  GOTO  720 

MAXTRW  = DTRW 

720  DTRW  = MAXTRW 

721  TRWDUM  = MAXTRW 
C 

YRMAX(IP,INM)  = DYMAX 
TRWR (IP, INM)  = DTRW 
TWR(IP,INM)  = TRWDUM 
MAXTRWR(IP, INM)  = MAXTRW 
SPCR  = DUMYF  * PC4 
DAMR  = BD  * ABD  * YRD  * PCS 
FDR  = SPCR  - DAMR 
C 

C CALCULATE  SPRING  FORCES 
C 


C 


DPHIS 

DPHISD 

COSPHIS 

SINPHIS 

CDPHIS 

SDPHIS 

TDPHIS 

CDPHIS2 


PHIS  - PHIU 
PHISD  - PHIUD 
DCOS (PHIS) 
DSIN(PHIS) 

DCOS (DPHIS) 
DSIN( DPHIS) 
DTAN( DPHIS) 
CDPHIS  * CDPHIS 


DELZ  = ZUP  - ZSP 
DELZD  = ZUDP  - ZSDP 
DELY  = YS  - YU 
DELYD  = YSDP  - YUDP 
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c 


DELSQ  = DSQRT(DELY*DELY  + DELZ*DELZ) 

DELSQD  = (DELZ*DELZD  + DELY*DELYD)  / DELSQ 
C 

C CALCULATE  LEFT  SPRING  AND  BUMP  STOP  FORCE 
C 

SLL  = (DELSQ+HU)  / CDPHIS  - HS  + S * TDPHIS 

SLLD  = DELSQD/CDPHIS  + DPHISD* (( DELSQ+HU) *SDPHIS+S ) /CDPHIS2 

DEFLL  = SPRLNG-SLL 

IF  (SLL.LE.BSSIZL)  GOTO  553 

IF  (SLL.GE.BSSIZU)  GOTO  555 

AK2  = O.ODO 

AB2  = O.ODO 

GOTO  558 

553  AK2  = l.ODO 

AB2  = l.ODO 
DBL  = BSSIZL  - SLL 
GOTO  556 

555  DBL  = BSSIZU  - SLL 

556  DBLD  = SLLD 
C 

558  FSL  = SUSKI  * DEFLL  - B1  * SLLD  + 

1 AK2  * SUSK2  * DBL  - AB2  * B2  * DBLD 

C 

C CALCULATE  RIGHT  SPRING  AND  BUMP  STOP  FORCE 
C 


572 


574 

575 
C 

577 

C 


C 

C 


SLR  = (DELSQ+HU)  / CDPHIS  - HS  - S * TDPHIS 

SLRD  = DELSQD/CDPHIS  + DPHISD* ( (DELSQ+HU) *SDPHIS-S) /CDPHIS2 

DEFLR  = SPRUNG  - SLR 

IF  (SLR. LE. BSSIZL)  GOTO  572 

IF  (SLR. GE. BSSIZU)  GOTO  574 

AK2  = O.ODO 

AB2  = O.ODO 

GOTO  577 

AK2  = l.ODO 

AB2  = l.ODO 

DBR  = BSSIZL  - SLR 

GOTO  575 

DBR  = BSSIZU  - SLR 
DBRD  = SLRD 

FSR  = SUSKI  * DEFLR  - B1  * SLRD  + 

1 AK2  * SUSK2  * DBR  - AB2  * B2  * DBRD 

SDPHIU  = DSIN(PHIU-PHIS) 

CDPHIU  = DCOS (PHIU-PHIS) 

DPHIUD  = PHIUD-PHISD 

DELP  = DELY*COSPHIS  - DELZ*SINPHIS  - (HR-HU) *SDPHIU 
DELPDl  = PHISD  * (DELZ*COSPHIS  + DELY*SINPHIS) 

DELPD2  = DELYD  * COSPHIS  - DELZD  * SINPHIS 
DELPD3  = (HR-HU)  * CDPHIU  * DPHIUD 
DELPD  = DELPD2  - DELPDl  - DELPD3 

PINFl  = DELP*SPRNGK  + BL*DELPD 
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CALCULATE  THE  INTEGRATION  ARGUMENTS 


YSDD  = ( (FSL+FSR) *SINPHIS  - PINFl*COSPHIS)/MS 

ZSDD  = (-(FSL+FSR) *COSPHIS  - PINF1*SINPHIS)/MS  + G 

HPl  = DELY  - HR*SINPHIU 
HP2  = DELZ  - HR*COSPHIU 
HPINFl  = DSQRT(HP1*HP1  + HP2*HP2) 


PHISDD  = ((FSL-FSR)  * S + PINFl  * HPINFl)  / IS 

FYL  = FYLF  + FYLR 
FYR  = FYRF  + FYRR 
FZL  = FZLF  + FZLR 
FZR  = FZRF  + FZRR 

YUDD  = (PINFl*COSPHIS- (FSL+FSR) +SINPHIS- (FYR+FYL+FDF+FDR) ) /MU 
ZUDD  = ( (FSL+FSR) *COSPHIS+PINFl*SINPHIS-FZL-FZR) /MU  + G 


1 
1 

PHIUDD  = (TERM1+TERM3-TERM4-TERM5-TERM6+TERM7+TERM9)/IU 
THEDD  = ( (IZ-IX) *YAWR*PHISD  - B* (FZLR+FZRR) +A* (FZLF+FZRF) ) /lY 


FREDLF  = (HU  + DABS(ZTLF))  * SINPHIU 
FREDRF  = (HU  + DABS(ZTRF))  * SINPHIU 
FREDLR  = (HU  + DABS(ZTLR))  * SINPHIU 
FREDRR  = (HU  + DABS(ZTRR))  * SINPHIU 
BOBl  = -ZUP  - HCURB 


TERMl 

TERMS 


= (FSR-FSL)  * S * DCOS(DPHIS) 


FZLF  * (TRW/2. ODO  * 
FZLR  * (TRW/2. ODO  * 
TERM4  = FZRF  * (TWF(IP,INM) 
FZRR  * (TWR(IP,INM) 
TERMS  = FYL  * ZUP 
TERM6  = FYR  * ZUP 
TERM?  = (HR  - HU)  * PINFl  * 
TERM9  = FDF  * BOBl  + FDR  * 


COSPHIU  + FREDLF) + 
COSPHIU  + FREDLR) 

* COSPHIU  - FREDRF) + 

* COSPHIU  - FREDRR) 


DCOS (DPHIS) 
BOBl 


YAWRD 

1 


((FYLR  + FYRR  + FDR) *B  - (FYLF  + FYRF  + FDF) *A  + 
PHISD  * THED  * (IX  - lY) ) / IZ 


YPRIME(l) 
YPRIME(2) 
YPRIME(3) 
YPRIME(4) 
YPRIME(5) 
YPRIME(6) 
YPRIME(7) 
YPRIME ( 8 ) 


= YSDD 
= ZSDD 
= PHISDD 
= YUDD 
= ZUDD 
= PHIUDD 
= THEDD 
= YAWRD 


YPRIME (9)  = YSD 

YPRIME (10)  = ZSD 


ooo  ooo  ooo  o ooo  ooo  ooo  ooo 


c 

c 

c 

c 


YPRIME(ll) 

YPRIME(12) 

YPRIME(13) 

YPRIME(14) 

YPRIME(15) 

YPRIME(16) 

YPRIME(17) 


PHISD 

YUD 

ZUD 

PHIUD 

THED 

YAWR 

YUD*COSPSI+YAWR* (A*COSPSI+ 
TRW/2. ODO*SINPSI) 


SET  VARIABLE  VALUES  EQUAL  TO  VALUES  CALCULATED  BY  SUBROUTINE 
MODEL 


VAR (18)  = V 
VAR (19)  = FZL 
VAR (20)  = YPRIME(3) 
VAR(21)  = YPRIME(6) 


CALCULATE  THE  MASS  MOMENTS  OF  INERTIA  AROUND  PT  A 


YUA  = (TRAD  + HU)  * SINPHIU  - 
1 (TRWF(IP,INM)  + TRWR(IP,INM) ) / 2 . ODO  * COSPHIU 

YSA  = YUA  + (YS  - YU) 

ISA  = IS  + MS  * ( (ZSP  + HCURB)  * (ZSP  + HCURB)  + YSA* YSA  ) 
lUA  = lU  + MU  * ( (ZUP  + HCURB)  * (ZUP  + HCURB)  + YUA*YUA  ) 

CALCULATE  THE  INSTANTANEOUS  AVG  C.G.  HEIGHT 

H3  = “(MS*ZSP  + MU*ZUP)  / (MS+MU) 

CALCULATE  THE  NONCENTRIODAL  K.E. 

TNR  = (ISA  * PHISD*PHISD  + lUA  * PHIUD*PHIUD)  / 2 . ODO 
VAR (22)  = TNR 

CALCULATE  THE  GENERAL  PLANE  MOTION  ENERGY 

TGPM  = ( MS* (YSDP*YSDP  + ZSDP*ZSDP)  + MU* ( YUDP*YUDP  + ZUDP*ZUDP) 
1 + IS*PHISD*PHISD  + IU*PHIUD*PHIUD  ) / 2 . ODO 

VAR (23)  = TGPM 

CALCULATE  THE  TRANSLATIONAL  ENERGY 

TRAN  = TGPM  - TNR 
VAR (24)  = TRAN 

CALCULATE  THE  POTENTIAL  ENERGY  OF  THE  VEHICLE 

VENG  = (MS  + MU)  * G * (H3  - H4(IP,INM)) 

VAR(25)  = VENG 

CALCULATE  THE  TOTAL  ENERGY  OF  THE  SYSTEM 

VAR (26)  = VENG  + TGPM 
C 
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CALCULATE  THE  CRITICAL  ENERGY  OF  ROLOVER 

SRSF  = TRW  / 2.0D0  / H4(IP,INM)  * PCI 
PCRIT  = G * (MS  + MU)  * H4(IP,INM)  * ( DSQRT(1.0D0 
1 - l.ODO  ) 

VCRIT  = PCRIT 

VAR (27)  = PCRIT  - TNR 

UPDATED  SEC.  PARA  SENT  TO  DGEAR 

CALL  VDATA 

REASSIGN  VARIABLE  VALUES 


Y(l) 

= 

YSD 

Y(2) 

= 

ZSD 

Y(3) 

= 

PHISD 

Y(4) 

= 

YUD 

Y(5) 

= 

ZUD 

Y(6) 

= 

PHIUD 

Y(7) 

= 

THED 

Y(8) 

= 

YAWR 

Y(9) 

= 

YS 

Y(10) 

= 

ZS 

Y(ll) 

= 

PHIS 

Y(12) 

= 

YU 

Y(13) 

= 

ZU 

Y(14) 

= 

PHIU 

Y(15) 

= 

THE 

Y(16) 

= 

PS  I 

Y(17) 

YGL 

RETURN 

END 

SRSF**2) 
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Data  File  for  Sensitivity  Program 
1975  Nova 


1, 

1, 

1, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

1, 

0, 

1, 

0, 

0, 

0, 

1, 

0, 

0, 

0, 

0, 

0, 

0, 

1, 

1, 

1, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 


DIST  FROM  C.G.  TO  FRONT  AXLE 
WHEELBASE 

VEHICLE  TRACK  WIDTH 

DIST  FROM  AXLE  TO  UNSPRUNG  MASS  CG 
SINGLE  TIRE  COMPRESSION  RATE  CONST. 
SINGLE  TIRE  DAMPING  CONST. 

COEFF.  OF  FRICTION  TIRE/ GROUND 
CURB  FORCE  DAMPING  CONSTANT 
SLOPE  OF  WHEEL  COMPLIANCE  AT  CURB 
SLOPE  OF  CAR  DEFORMATION  AT  CURB 
TOTAL  UNSPRUNG  MASS 

DIST  FROM  AXLE  TO  LAT.  REST.  FORCE 
SUSPENSION  SPRING  CONST 
BUMP  STOP  SPRING  CONST 
SUSPENSION  DAMPING  CONST. 

BUMP  STOP  DAMPING  CONST. 

TOTAL  ROLL  INERTIA 
TOTAL  PITCH  INERTIA 
TOTAL  YAW  INERTIA 

DIST  FROM  SPR.  MASS  BOTT.  TO  SMCG 
TOTAL  UNSPRUNG  MASS  ROLL  INERTIA 
CONSTANT  FORCE  OF  WHEEL  DEF.  0 CURB 
VEL.  OF  CONSTANT  SIDE  VEL.  FUNCTION 
CRUSH  DIST  WHERE  SUSPENSION  DEFORMS 
CRUSH  DIST  WHERE  BODY  DEFORMS 
SPRUNG  MASS 

SPRUNG  MASS  ROLL  INERTIA 
TIRE  RADIUS 
CURB  HEIGHT 

HALF  DIST  OF  SPRING  BASE 
LATERAL  BUSHING  SPRING  CONST. 
UNLOADED  SUSPENSION  SPRING  LENGTH 
DIST  SS  AXLE  TO  BUMP  STOPS 
LATERAL  BUSHING  DAMPING  CONSTANT 
LATERAL  POSITION  OF  CURB 
DUMMY  TEST  PARAMETER 
ENDPP 


, 4.45D0 
,9.25D0 
,5.0D0 
,0.2D0 
, 12000. ODO 
,80.0D0 
, 0.75D0 
, 0.25D0 
,24000. ODO 
,96000. ODO 
, 16.7D0 
, l.ODO 
,6900. ODO 
, 33000. ODO 
, 325. ODO 
,33. ODO 
,330. ODO 
, 2240. ODO 
,2240. ODO 
,1.2D0 
, 132 . ODO 
,4000. ODO 
, 1. ODO 
, 0. 1667D0 
,0.6667D0 
, 105. ODO 
,240. ODO 
, l.ODO 
,0.5D0 
, 1.75D0 
,200000. ODO 
,0.6667D0 
, 0.25D0 
,1000. ODO 
,7.50D0 
, 1. ODO 
, 0. ODO 
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Data  File  for  Sensitivity  Programs  (Con't) 

1975  Nova 


0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

1/ 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 


PCI  - SRSF 

' , 1. ODO 

PC2  - SPCF 

' , 1. ODO 

PC 3 - DAMF 

' , 1. ODO 

PC4  - SPCR 

' , 1. ODO 

PC5  - DAMR 

' , 1. ODO 

PC6  - ZSP 

' , 1. ODO 

PC7  - ZUP 

' , 1. ODO 

PCS  - ZSDP 

' , 1. ODO 

PC9  - ZUDP 

' , 1. ODO 

PCIO  - YSDP 

' , 1. ODO 

PC 11  - YUDP 

' , 1. ODO 

ENDSP 

' , 0. ODO 

1 YS  Dot 

’ , O.ODO 

2 ZS  Dot 

' , 0. ODO 

3 PHIS  Dot 

' , O.ODO 

4 YU  Dot 

' , 0. ODO 

5 ZU  Dot 

' , O.ODO 

6 PHIU  Dot 

' , 0. ODO 

7 THETA  Dot 

' , 0. ODO 

8 r 

' , 0. ODO 

9 YS 

' , 0. ODO 

10  ZS 

' , 0. ODO 

11  PHIS 

' , 0. ODO 

12  YU 

' , 0. ODO 

13  ZU 

' , 0. ODO 

14  PHIU 

' , O.ODO 

15  THETA 

' , 0. ODO 

16  PSI 

' , 0. ODO 

17  Y GLOBAL 

' , 0. ODO 

18  V 

' , 0. ODO 

19  FZL 

' , 0. ODO 

20  PHISDD 

' , 0. ODO 

21  PHIUDD 

' , 0. ODO 

22  TNR 

' , 0. ODO 

23  TGPM 

' , O.ODO 

24  TRAN 

' , 0. ODO 

25  VENG 

' , 0. ODO 

26  TOTAL  ENERGY 

' , 0. ODO 

27  RPER 

' , 0. ODO 

ENDOV 

' , O.ODO 
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APPENDIX  E 


Parameter  Set  Description 
Geometric  Parameter  Set  #1: 


Label 

Description 

TRW 

hu 

Vehicle  track  width 

Height  of  unsprung  mass  c.g.  position  above  axle 
centerline 

hr 

Height  of  lateral  restraining  force  member  (roll 
axis)  above  the  axle  centerline  (average  front  and 

hs 

rear) 

Static  height  of  sprung  mass  c.g.  position  above 
axle  centerline 

L 

a 

Vehicle  wheelbase 

Horizontal  distance  from  the  vehicle  c.g.  position 
to  the  front  axle 

SRSF 

The  ratio  of  track  width  to  c.g.  height 

Geometric  Parameter  Set  #2: 


Label 

Description 

^curb 

SPRLNG 

T 

ABSDEF 

Height  of  the  curb 

The  undeflected  length  of  suspension  springs 
Undeformed  tire  radius 

Vertical  wheel  travel  distance  from  the  static 

s 

position  of  the  wheel  to  the  point  where  contact 
with  the  suspension  bump  stop  begins 
Half  of  the  distance  between  left  and  right  suspen- 
sion springs 

Yu 

Initial  distance  between  vehicle  and  curb  (measured 
from  front  inside  tire) 

E-1 


Mass  Parameter  Set: 


Label 

Descrintion 

M 

^^s 

Mu 

Ix 

Total  sprung  mass 
Total  unsprung  mass 

Total  roll  inertia  about  total  mass  center  of 
gravity 

Iy 

iz 

Is 

lu 

Total  pitch  inertia 
Total  yaw  inertia 

Sprung  mass  roll  inertia  about  its  center  of  gravity 
Total  unsprung  mass  roll  inertia  about  center  of 

gravity 

Stiffness  Parameter  Set: 


Label 

Descriotion 

Suspension  spring  rate  per  side,  equivalent  acting 
on  wheel 

^2 

Bump  stop  springrate  per  side,  equivalent  acting  at 
wheel 

K 

Kz 

KDl 

Bushing  (lateral  restraining  force)  spring  constant 
Single  tire  radial  spring  rate 

Wheel  and  suspension  lateral  deflection  spring  rate 
when  striking  the  curb 

KD2 

Slope  of  force  vs.  lateral  deformation  of  body, 
frame  and  engine  when  striking  the  curb 

Damping  Parameter  Set: 


Label 

Descriotion 

Suspension  shock  absorber  damping  rate  per  side, 
equivalent  acting  on  wheel 

B2 

Bump  stop  damping  rate  per  side,  equivalent  acting 
at  wheel 

B 

Bd 

Bushing  (lateral  restraining  force)  damping  constant 
Single  tire  radial  damping  rate 

Vehicle  structure  damping  constant  in  the  lateral 
deformation  direction 
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Force  Parameter  Set: 


Label 

WHLDEF 

Kv 

POINTl 

P0INT2 


Description 


Coefficient  of  friction  between  tire/ground 
Constant  force  during  crush  deformation  of  the 
"soft"  suspension  and  wheel  members  when  striking  a 
curb 

Velocity  at  which  tire  slides  with  a constant  side 
velocity  function 

Deflection  distance  at  which  wheel  and  suspension 
begin  to  crush  (permanent  deformation) 

Crush  distance  at  which  the  more  solid  members  of 
the  body,  frame,  and  engine  become  involved  in  lateral 
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